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namic processes in systems of interacting particles is formulated on the basis 
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the Enskog-Landau kinetic equation for a system of charged hard spheres. The 
BBGKY hierarchy is analyzed on the basis of modified group expansions. Gen- 
eralized transport equations are obtained in view of a self-consistent description 
of kinetics and hydrodynamics. Time correlation functions, spectra of collective 
excitations and generalized transport coefficients are investigated in the case of 
weakly nonequilibrium systems of interacting particles. 

Key words: kinetics, hydrodynamics, kinetic equations, transport coefficients, 
(time) correlation functions 

PACS: 05.20.Dd, 05.60. +w, 52.25.Fi, 71.45.G, 82.20.M 



1. Introduction 

The problem of a self-consistent description of fast and slow processes which are 
connected with both linear and non-linear fluctuations of observed quantities for 
various physical systems, such as dense gases, liquids, their mixtures and plasma, 
remains an actual problem in nonequilibrium statistical mechanics. An important 
task in this direction is the construction of kinetic equations for dense systems with 
taking into account collective effects (hydrodynamic contribution) into the collision 
integrals. On the other hand, there is a problem of the calculation of time cor- 
relation functions and collective excitations spectra in the range of intermediate 
values of wavevectors and frequencies, since low- and high-frequencies regions are 
described sufficiently well within the frameworks of molecular hydrodynamics and 
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kinetic equations, respectively. In the intermediate region, the kinetic and hydrody- 
namic processes are connected and must be considered simultaneously. 

A significant interest in these problems was exhibited by Prof. D.N.Zubarev in 
his investigations. One of the approaches to the unification of kinetics and hydro- 
dynamics in the theory of transport phenomena for systems of interacting particles 
was proposed by Zubarev and Morozov in [ |IJ and developed in further papers [ 
H, H, H, H, H, 0, |, H, |l(J [□], |l2i . This statistical approach is based on a modifica- 
tion of the boundary conditions for the weakening of correlations to the BBGKY 
hierarchy for many-particle distribution functions and to the Liouville equation for 
a full nonequilibrium distribution function (nonequilibrium statistical operator - 
NSO). The modification of the boundary conditions consists in the fact that both 
the one-particle nonequilibrium distribution function and the local conservation laws 
for mass, momentum and total energy are included into the parameters for an ab- 
breviated description of a nonequilibrium state of the system. The concept of a 
self-consistent description of kinetics and hydrodynamics was applied to plasma in 
an electro-magnetic field [ [13], [14], [T5|], to quantum systems with taking into account 
coupled states [ ||, [10], [T 
thermo field dynamics [ [l 



quantum Bose systems [ [17]] and the nonequilibrium 
The questions as to the necessity and possibility of a 
unified description of kinetic and hydrodynamic processes were discussed in papers 



by Klimontovich [ |lg |g, |l 



In this paper we present some actual results obtained recently within the concept 
of a self-consistent description of kinetics and hydrodynamics for dense gases and 
simple liquids. In section 2 we present the results for kinetic equations of RET [ 
|3], 0, GDRS [ §] and Enskog-Landau theories [ 0, 0], obtained in the polarization 
approximation for a system of hard spheres and a Lenard-Balescu kinetic equation 
[ [TIL |12| derived from the BBGKY hierarchy with modified boundary conditions. 
Generalized transfer equations, time correlation functions and transport coefficients 
are presented in section 3 on the basis of a self-consistent description of kinetics 
and hydrodynamics within the framework of the nonequilibrium statistical operator 
method. 



2. Conception of a consistent description of kinetics and hydro- 
dynamics for dense gases and liquids 

2.1. Overview 

Bogolubov [ 0] proposed a consequent approach to the construction of kinetic 
equations which is based on the chain of equations for s-particle distribution func- 



tions and a weakening correlation principle. There are a lot of theories [|23], ^5], 
which differ in forms of presentation, but all these approaches use the same weaken- 
ing correlation condition and are most effective when a small parameter (interaction, 
density, etc.) is presented. 

In the kinetic theory of classical gases we can point out two principal problems. 
The first problem is connected with the fact that collision integrals can depend 
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nonanalytically on density. Because of this, even in the case of a small-density gas, 
in order to calculate corrections to transport coefficients it is necessary to perform 
a partial resummation of the BBGKY hierarchy [ 27, 28fl . The second problem is 



the construction of kinetic equations for the gases of high densities. In this case 
we cannot restrict ourselves by several terms of the expansion of collision integrals 
on densities and the analysis of the BBGKY hierarchy becomes very complicated 
(generally speaking, the density is not a small parameter here). That is why the 
construction of kinetic equations for dense gases and liquids with model interparticle 
potentials of interaction has a great importance in the problem under investigation. 

The first theory in this direction is a semi-phenomenological standard Enskog 
theory (SET) [ |2^, Q of dense gases. Ideas similar to those used at the derivation 
of the Boltzmann equation were also used in this approach. The Enskog equation 
plays an essential role in the kinetic theory, next to the Boltzmann one [ |30| . This 
equation was obtained at the modelling of molecules by hard spheres for dense 
gases. As a result, the collision integral was presented in an analytical form. It 
is obtained by means of a hard sphere model where collisions can be considered 
as momentary and by the fact that the multiparticle contact at the same time is 
reputed to be infinitely small. Density correction introduced by Enskog proved to 
be considerable, as far as transport, due to collisions in a dense system, is the main 
mechanism of transport. Each molecule is almost localized at one point of space by 
the surrounding neighbour molecules and, therefore, the flow transport is suppressed. 
Though this theory properly describes density dependence of the kinetic coefficients, 
at the same time the suppositions about the structure of the collision integral remain 
sufficiently rough and phenomenological. Despite approximate assumptions on the 
collision integral in the kinetic equation for a one-particle distribution function of 
hard spheres, the Enskog theory very well describes a set of properties for real dense 
gases [H, gj. 

Davis, Rise and Sengers [ [32[] proposed a kinetic DRS theory. Within the frame- 
work of this theory, the interparticle potential of an interaction is chosen in a square 
well form. An attractive part of the real potential is approximated here by a finite 
height wall. 

From the point of view of the statistical theory of nonequilibrium processes, 
SET and DRS contain two essential drawbacks. The first one is that their kinetic 
equations are not obtained within the framework of some consequential theoretical 
scheme and one does not know how to improve these theories. And the second one 
is that the if-theorem has not been proved. Nevertheless, not long ago a way of 
constructing the SET entropy functional was given [ In order to overcome these 
drawbacks the authors of [ |34| obtained a kinetic equation for a revised version of 
the Enskog theory (RET) using the diagram method. Resibois [ proved the 

//-theorem for it. In 1985, a revised version of the DRS theory (RDRS) was proposed 
37j . The kinetic equation of RDRS satisfies the if-theorem as well. A generalized 



version of RDRS — GDRS was considered in [ @, , being based on the approaches 
of [ HJ, 0. Here, in order to treat a more realistic model, a kinetic equation for dense 
classical systems is offered with an interparticle interaction potential in the shape 
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of a multistep function, where its sequential derivation and the normal solution are 
obtained. 



In the paper by Stell et al. [ |38], |39j the authors completed the construction of 
the kinetic variational theory (KVT) introducing a whole family of theories, namely, 
KVT I, KVT II and KVT III. The KVT III refers to the local energy constraint 
in maximizing entropy as distinguished from the global energy constraint (KVT 
II) or the hard-core constraint (KVT I). According to this classification, the kinetic 
equation of RDRS [ [37] is obtained by applying KVT III to the square well potential 
and KMFT [ [3£| is derived from KVT I as its application to the potential with a 
smooth tail. In [ [591, the KVT III version of KMFT was proposed. In this case 
the quasiequilibrium binary distribution function (QBDF) of hard spheres is to be 
substituted by the full QBDF, which takes into account the full potential of the 
interaction. The main conclusion of the KVT I version of KMFT [ f40| is that the 
smooth part of the potential does not contribute explicitly to transport coefficients. 
This fact is caused by the approximations of the theory. When the full QBDF is used, 
as it is in the KVT III version of KMFT, transport coefficients are determined by the 
soft part of the potential. The case when the potential between the hard sphere and 
the attractive walls is considered as a smooth tail (instead of the constant potential 



in DRS) has also been investigated [ 41 1. The kinetic equation was obtained for this 
potential applying a mean field approximation for a smooth tail. 

At this stage of investigations in the kinetic theory of dense gases and liquids, 
paper [ |T| by Zubarev and Morozov played a fundamental role. In this paper a new 
formulation of boundary conditions for the BBGKY hierarchy is presented. Such a 
formulation took into account correlations connected with local conservation laws. 
In the binary collision approximation, this modification of the weakening correlation 
condition by Bogolubov led to a kinetic equation for a system of hard spheres which 
is close in structure to the usual Enskog kinetic equation. Similar ideas were proposed 
independently by the authors of paper [ 37 1 at the derivation of a kinetic equation 
when an interaction between the particles is modelled by a square well potential. It 
is necessary to point out that a somewhat different modification of the Bogolubov 
approach [ |22| was presented in the papers by Rudyak [ |42], |3[ [0J] . This allows one 
to obtain an Enskog-like kinetic equation and extend it to the systems with a soft 
potential of interaction between the particles. 

The modification of weakening correlation conditions to a chain of equations by 
Bogolubov was developed in papers [ ||, |3|, f|, [5], || [7|, §, |9j, [TTJ, [12] . An important 
achievement of the given approach is the fact that the modified boundary condition 
gives the possibility [ f|, |7|, [L4|] to derive consequently the kinetic equation of the 
RET theory [ |34| for the first time. On the basis of result [|lj, a kinetic equation was 
obtained [ 0, || for systems with a multistep potential of interaction (in particular, 
the if -theorem was proved for this equation in paper [ |J), whereas an Enskog- 
Landau kinetic equation was derived [ [3|, [7| for a system of charged hard spheres. 

Normal solutions for the obtained kinetic equations were found using the Chap- 
man-Enskog method. With the help of these solutions, the numerical calculations 
of transport coefficients (bulk and shear viscosities, thermal conductivity) were car- 
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ried out for Argon-like systems [ [2|, 45fl , ionized Argon [ [7| and mixtures [ 46|. It is 
important to stress that the RET kinetic equation, the kinetic equation for systems 
with multistep potentials of interaction and the Enskog-Landau kinetic equation 
for charged hard spheres were derived from the BBGKY hierarchy with a modified 
boundary condition in a binary collision approximation. Obviously, the presented 
equations have a restricted region of application. They cannot be applied to the de- 
scription of systems with significantly collective effects caused by Coulomb, dipole 
or other long-range forces of interaction between the particles. To describe the col- 
lective effects in systems with a long-range character of interaction it is necessary to 
consider higher-order approximations for collision integrals. An example of such an 
equation is the Lenard-Balescu kinetic equation [ ^7\, |48], for Coulomb plasma. 



To analyze solutions to the BBGKY hierarchy with a modified boundary condi- 
tion in higher-order approximations on interparticle correlations we have applied a 



concept of group expansions [ [U], [12 



2.2. The Liouville equation and the BBGKY hierarchy with a modified 
boundary condition 

Let us consider a system of N identical classical particles which are enclosed in 
volume V, with the Hamiltonian: 

N 2 -, N N 



2m 2 

3=1 j=l k=l 

where $ (|fjjfe|) is the interaction energy between two particles j and k; \rjk\ = \rj~ 
is the distance between a pair of interacting particles; Pj is the momentum of jth 
particle and m denotes its mass. A nonequilibrium state of such a system is described 
by the A-particle nonequilibrium distribution function g (x N ; t) = q (xi, . . . , x^; t) 
which satisfies the Liouville equation: 

^+iL N ^g(x N ,t) =0, (2.2) 

where i = ^/—l, Xj = {rj,Pj} is a set of phase variables (coordinates and momenta), 
L N is the Liouville operator: 

N N N 

j=l j=l k=l 
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where 

m dr -i 



L(j,k)= i 



M{\r jk \) ( d d 



dr jk \dpj dp k/ 

The function g {x N \ t) is symmetrical with respect to permutations x\ ^ Xj of 
phase variables for an arbitrary pair of particles. This function satisfies the normal- 
ization condition 

J dr N g (x N ; t) = 1, dr N = (dx) N /N\ = drdp. (2.4) 



In order to solve the Liouville equation Q2.2|), it is necessary to introduce a 
boundary condition. It must be chosen in such a way that the solution to the Liouville 
equation will correspond to a physical state of the system under consideration. Using 
the nonequilibrium statistical operator method by D.N.Zubarev pi], [5C], |51| (NSO 
method), we shall search for such solutions to the Liouville equation depending on 
time explicitly via the values for some set of observable variables, which is sufficient 
for describing a nonequilibrium state without depending on the initial time t Q . 

The solution to the Liouville equation, which satisfies the initial condition 

g (x N ;t) \ t=to = g q (x N ;t ) , 

has the form: 

g(x N ;t,t )=e- iL ^-^g q (x N ;t ). (2.5) 

We shall consider such times t ^> t when the details of the initial state of the 
system can be neglected. Then, to avoid the dependence on t , let us average ( |2.5|) 
with respect to the initial times from to t t and perform the boundary transition 
(t — t Q ) — >• oo. As a result, one obtains [ |9], [10], |49| : 



g(x N ;t)=e J dt' e £t 'e iLNt 'g q (x N ;t + t') , t' = t - t, 



(2.6) 



where e tends to +0 after the thermodynamic limit transition. It can be shown by 
straightforward differentiations that solution ( |2.6| ) satisfies the Liouville equation 
with an infinitesimal source in the right-hand side: 



^ + iL N ^j g (x N , t) = -e (q (x N , t) - g q (x N , t) 



(2.7) 



The source breaks the symmetry of the Liouville equation with respect to time inver- 
sion and selects retarded solutions which correspond to an abbreviated description 
of a nonequilibrium state of the system. The auxiliary function £ q (x N ;t) denotes 
a quasiequilibrium distribution function that is defined from an extreme condition 
for the informational entropy of the system, provided the normalization condition 
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is preserved and the average values for variables of the abbreviated description are 
fixed. 

The choice of g q (x N ] i) depends mainly on a nonequilibrium state of the system 
under consideration. In the case of low density gases, where times of free motion are 
essentially larger than collision times, higher-order distribution functions of particles 
become dependent on time only via one-particle distribution functions [ [E| |52|, |53| . 
It does mean that an abbreviated description of nonequilibrium states is available, 
and the total nonequilibrium distribution function depends on time via fi(x;t). In 
such a case, the quasiequilibrium distribution function g q (x N ;t) reads [|52|, |53|j : 

g,(x N ;t)=f[^^, (2.8) 

3=1 6 

where e is the natural logarithm base. Then, the Liouville equation with a small 
source (|2.7|) in view of (|2.8| ) corresponds to the abbreviated description of time 
evolution of the system on a kinetic stage, when only a one-particle distribution 
function is considered as a slow variable. However, there are always additional quan- 
tities which vary in time slowly because they are locally conserved. In the case of 
a one-component system, the mass density p(r;t), momentum j(r;t) and total en- 
ergy S(r;t) belong to such quantities. At long times they satisfy the generalized 
hydrodynamics equations. Generally speaking, the equation for fi(x;t) must be 
conjugated with these equations. For low density gases, such a conjugation can be 
done, in principle, with an arbitrary precision in each order on density. In high 
density gases and liquids, when a small parameter is absent, the correlation times 
corresponding to hydrodynamic quantities become commensurable with the char- 
acteristic times of varying one-particle distribution functions. Therefore, in dense 
gases and liquids, the kinetics and hydrodynamics are closely connected and should 
be considered simultaneously. That is why many-particle correlations, related to the 
local conservation laws of mass, momentum and total energy, cannot be neglected 
[ EO, H, H llOU - The local conservation laws affect the kinetic processes due to an 
interaction of the selected particles group with other particles of the system. This 
interaction is especially important in the case of high densities, and it must be taken 
into consideration. Then, at the construction of kinetic equations for high densities, 
it is necessary to choose the abbreviated description of a nonequilibrium system in 
the form to satisfy the true dynamics of conserved quantities automatically. To this 
end, the densities of hydrodynamic variables should be included together with the 
one-particle distribution function fi (x; t) into the initial set of parameters of the 
abbreviated description [ 0]j @> [T0[1 - The next phase functions correspond to the 
densities of the hydrodynamic variables p(r; t), j(r; t) and S(r; t): 



7 



M.V.Tokarchuk, I.P.Omelyan, A.E.Kobryn 



p(r) = J dp ni(x)m, 

j(r) = j dp ni(x)p, (2.9) 

f P 2 1 f 

Sir) = / dpnAx) 1 — / dr dp dp' h 2 (x, x')<P(\r — r I), 

J 2m 2 J 

where h\(x) and n2(x,x') are one- and two-particle microscopic phase densities by 
Klimontovich [|||: 

N N 

ni(x) =^26(x-x j ) = ^25(r-r j )8(p-p j ), (2.10) 

3=1 3=1 
N 

5{x-Xj)5{x' -x k ). (2.11) 



n 2 (x, x' 



Relations (|2.9|) show a distinctive role of the potential interaction energy. Contrary to 
p(r; t) = (p(r)Y and j(r; t) = (j(r)Y, the nonequilibrium values of the total energy 
£(r;t) = (S(r)Y cannot be expressed via the one-particle distribution function 
fi(x] t) = (fii(x)Y only, because in order to evaluate the potential part of S int (r; t) = 
(£mt( r )Y it is necessary to involve the two-particle distribution function f2(x, x'\ t) = 
(77-2(0;, x')Y- Here 

4t(r) = \J dr ' d P V Mx, x')$(\r - r'|) (2.12) 

is the density of the potential energy of interaction. The next conclusion can be 
formulated as follows. If the one-particle distribution function fi(x; t) is chosen as a 
parameter of an abbreviated description, then the density of the interaction energy 
( |2.12| ) can be considered as an additional independent parameter. One can find the 
quasiequilibrium distribution function g q (x N ; t) from the condition of extremum for 
a functional of the informational entropy S in f (t) = — f dr N g (x N ; t) In g (x N ; t) at 

fixed average values of (ni(x)Y = fi(x;t), (£ int (r))' = £ int (r;t), including the nor- 
malization condition for g q (x N ;t) [|^, [| |TU], f dr N g q (x N ;t) = f dr N g (x N ;t) 
= 1. This is equivalent to finding an unconditional extreme for the functional 



L{q) = J dr N g(x N ;t) ^l-1ng(x N ;t) - -J dr/3(r;t)£ mt (r) -J dxa(x;t)hi(x) 

where @(r;t), a(x;t) are the Lagrange multipliers. 

Taking the variation —L(g), after some simple manipulations one obtains the 

og 

quasiequilibrium distribution function 

g q (x N -i) = exp -J dr f3(r;t)S int (r) -J dx a{x;t)hi{x)^ , (2.13) 

$(t) = In / dr N expi-/ dr (3(r;t)S- mt (r) - dx a(x; t)fn(x) i. (2.14) 
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Here, is the Massieu-Planck functional. It is determined from the condition of 
normalization for the distribution g q (x N ;t). Relation fl2.13| ) was obtained for the 



first time in [[!]]. To determine the physical meaning of parameters /3(r; t) and a(x; t) 
let us rewrite g q (x N ; t) ( |2.13| ) in the form: 



g q (x N ;t) = expj-$(i) -J dr P(r;t)S'(r) - J dx a'(x; tfn^x) j , (2.15) 
= In J dr N exp j-^ dr 0(r;t)£'(r) - J dx a'(x; t)rii(x) 



where S'(r) is the density of the total energy in a reference frame which moves 
together with a system element of the mass velocity V(r; t) [ [5| 

S'(r) = S(r) - V(r; t)j(r) + -V 2 (r; t)n(r). (2.16) 

Here n(r) = J dp ni(x) is the density of the particles number. Parameters /3(r;t) 
and a'(x; t) in (|2.15|) are determined from the conditions of self-consistency, namely, 
the equality of quasi- average values (711(2;))* and (£'(r)) q with their real averages 

(nxCar))*, (£'(r))*: 

(£'( r ))t = {£'(r)Y, (...) q = J dr N ...g q (x ;t) . (2.17) 

In these transformations the parameters a'(x; t) and a(x; t) are connected by the 
relation 

a'(x; t) = a(x; t) - /3(r; f ) { £ - V(r; t)p + jV 2 (r; t) \ . 



In the case, when the conditions (|2.17|) take place, one can obtain some relations 



taking into account self-consistency conditions and varying the modified Massieu- 
Planck functional $(£) after ( 2.15|) with respect to parameters /3{r;t) and a'(x;t) 



-(£'(r)Y Q = -(E'(r)Y 



= -<ai(*))« = -<ni(x)>* = -h(x;t). 

It means that parameter /3(r; t) is conjugated to the average energy in an accompa- 
nying reference frame, and a'(x; t) is conjugated to the nonequilibrium one-particle 
distribution function fi(x;t). To determine the physical meaning of these parame- 
ters let us define the entropy of a system taking into account the self-consistency 
conditions ( |2.17|) : 



S(t) = -(\ng q (x N -t)Y q = <5>(t) + / dr /3(r; t)(S'(r)Y + / dx a' {x ] t){n l {x)) t . 

(2.19) 
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Taking functional derivatives of S(t) (|2.19|) with respect to (£'(r))' and (ni(x)Y at 
fixed corresponding averaged values gives the following thermodynamic relations: 

SS(t) P(r;t), WfL^affrt). (2.20) 



S(S'(r)Y v ' " 5fi(x;t) 

Hence, f3(r; t) is an analogue of the inverse local temperature. 

Two limiting cases follow from the structure of expression ( [2.19| ) for entropy. 
At a'(x;t) = —f3(r;t)n(r;t), ( |2.19| ) transforms into an expression for entropy which 
corresponds to the hydrodynamic description of the nonequilibrium state of the 
system [ §, g |TD[: 



S{t)=*{t) + J dr (3(r;t) ^(S'(r)Y - fi(r;t)(n(r)Y) , (2.21) 
with the following quasiequilibrium distribution function [ j5[ |9|, [10|1 : 

g q (x N ; t) = exp - J dr /3(r; t) (i\r) - ^(r; t)n(r)) } , (2.22) 

$(t) = In J dr N ex p|y dr P(r;t)(i'{r) - fi(r;t)h(r) 



and the corresponding self-consistency conditions for the definition of thermody- 
namic parameters /3(r; t), /i(r; t) (local value of the chemical potential): 

(n(r)>* q = (n(r)Y, (2'(r% = (S'(r))*. 

In the case when the contribution of the interaction energy between the particles 
can be neglected (a dilute gas), the quasiequilibrium distribution function ( [2.1 5|) has 
the form: 



g q (x N ;t) = exp |-$(t) - J dr (3{r; t)£( in (r) - J dx a'^t)^ 



x 



or, taking into account ( |2.16|) and the relation between a'(x;t) and a(x;t), one 
obtains: 

g q (x N ;t) =exp|-$(t)-y dx a(x;t)ni(x)\ , (2.23) 
$(i) = In J di~jv exp|— J dx a(x; t)fi\{x) 

where 

a a a 771 

%J?) = S kin (r) ~ V(r; t)j(r) + -V 2 (r; t)h(r) 



£kin(y) = J dp %-h\(x) is the density of the kinetic energy. Now, determining in 
( [2.23D parameter a(x;t) with the help of the self-consistency condition (n(x))* = 
(h(x)Y, one can show [ [5| that ( |2.23| ) for g q (x N ;t) transforms into distribution 
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Q2.8D when it is assumed that the only parameter of the abbreviated description for 
a nonequilibrium state of the system is a one-particle distribution function. As it 
is known [ [51], |53| , the quasiequilibrium distribution function (|2.23|) corresponds to 
the Boltzmann entropy of a dilute gas: 



s B (t) 



dx fi(x; t) In 



(2.24) 



In a general case, when kinetic and hydrodynamic processes are considered simul- 
taneously, the quasiequilibrium distribution function Q2.15|) or ( [2.131 ) can be rewrit- 
ten in a somewhat different form. This form is more convenient for the comparison 
with g q (x N ;t) ( |2.8|) , obtained in a usual way [ [22], |52|| , when fi(x;t) is the only 
parameter of the abbreviated description. First of all, let us note that one can in- 
clude parameter $(t) from ( p. 13 ) into parameter a(x; t) as a term which does not 
depend on x. Parameter a(x; t) in g q (x N ; t) can be excluded with the help of the self- 
consistency condition (h(x)) q = (n(x))* = fi(x;t). Reduction of g q (x N ',tj results 
in 



N 



0q (x N ;t)=exp{-U N (r N ;t)}H 



where functions u{rj\t) are obtained from the relations 

/j N—l N 
^-^e W {-U N (r,r^;t)}ll 



(2.25) 



u r,- 



3=2 



(2.26) 



1 N 

U N (r N ; t) = U N (n, ...,r N ;t) = - &(\ r i ~ r k\)P{r k ] t), 

n(r;t) = (n(r))* = j dp fi(x;t) is the nonequilibrium particles concentration. In 
expression (|2.25| ), C/jy(r; t) and u(r j] t) respectively, depend explicitly and implicitly 
on n(r;t) and /3(r;t) (or (S'(r)Y). To obtain the ordinary Bogolubov scheme [ 
22], |52]| , it is necessary to put Z7jy(r; t) = in ( |2.25| ) and (|2.26|) . Then, one can define 
it = e, and Q2.25 ) transforms into the quasiequilibrium distribution ( |2.8| ), as it should 
be. In a general case, u(r; t) is a functional of the nonequilibrium density of particles 
number n(r;t) and /3(r;t), which is an analogue of the inverse local temperature. 



Nevertheless, one should handle this analogy with far as definition ( |2.25|) 

can describe states which are far from local equilibrium. In particular, fi (x; t) can 
differ considerably from the local Maxwellian distribution. 

The entropy expression ( p,19| ) can be transformed according to the structure of 



the quasiequilibrium distribution function ( 2.25 ) 



S(t)= / dr P(r-t)(S int (r)Y 



dx fi (x\ t) In 



AM 

u(r;t) ' 



(2.27) 



Here the potential and kinetic parts are separated. In the case of low density gases, 
the influence of the potential energy can be neglected and u(r; t) = e. Then, expres- 
sion ( |2.27| ) tends to the usual Boltzmann entropy. 
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Thus, determining the quasiequilibrium distribution function g q (x N ;t) ( |2.25|) 



and entropy S(t) ( |2.27 ) of the system, when the nonequilibrium one-particle distri- 



bution function, as well as the average values of densities for the number of particles, 
momentum and energy are parameters of an abbreviated description of the nonequi- 
librium state which are locally conserved, the Liouville equation with the source 
can be presented in the form [ [l|. || : 



| + Q t)=s( e (x N ; t) - exp {-U N (r"; t) } ft 



(2.28) 

Further, on the basis of this equation, one obtains the BBGKY hierarchy for nonequi- 
librium distribution functions of classical particles with modified boundary condi- 
tions which take into account the nonequilibriumnes of the one-particle distribution 
function, as well as the local conservation laws. To obtain the first equation of the 
hierarchy, let us integrate with respect to the variables x N ~ x = {22, • • • , %n} the 
both parts of ( p.2S| ). Taking into account ( p.lOj ), ( |2.11|) and (|2.17|) , one obtains: 



0_ 

Of 



+ LL(1)^ fi{x!]t) + J dx 2 iL(l,2)/ 2 (zi,z 2 ;*)=0. 



(2.29) 



Integrating now ( 2.28| ) over the variables x N 2 = {x%, . . . ,x^}, after simple trans- 



formations one can obtain an equation for the two-particle distribution function 
f 2 (xx, X2] t) which differs from the corresponding equation in the Bogolubov hierar- 
chy [ [L], |3|] by a source in the right-hand side: 

(J^ + iL^ f 2 (x 1 ,x 2 ,t) + J djc 3 (iL(l,3) +iL(2,3) S jf 3 (x 1 ,X2,x 3 ;t) = 
-e(j 2 {xii x 2 ; t) - g 2 (n, r 2 ; t)/i(a?i; *)/i(2C 2 ; t)\ (2.30) 

In this equation, L 2 = L(l) + L(2) + L(l, 2) is the Liouville two-particle oper- 
ator and g2(ri,r2]t) denotes the binary coordinate distribution function for the 
quasiequilibrium state (|2.25|) ; f 2 (xx, x 2 ; t) is the nonequilibrium two-particle distri- 



bution function: 

f 2 (xi,x 2 ;t) = (n 2 (x 1 ,x 2 ;t)) t = / dr N _ 2 Q{xi,x 2 ,x N ~ 2 ;t). 



Similarly, integrating equation ( |2.28|) over x N ~ s = {x s+ \, . . . ,xn} phase variables, 
one derives next equations of the chain with the corresponding sources which define 
boundary conditions for reduced distribution functions. For the s-particle nonequi- 
librium distribution function 

f s (x s ; t) = (fisix 3 )) 1 = / dr N _ s q(x u ...,x s , x n ~ s ; t) 
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one has the following equation: 

^ + iL^/ s (x s ;t) + J dx s+1 J2iL(j,s + l)f s+1 (x s+1 ] t) = 

s 

-e(f s (x s ; t) - g s (r s ; t) J] f x { Xj \ tj) , (2.31) 



where 



s 1 s s 

L s = Y,L{j) + - L U,k), g s (r s ;t) = f s (r s ; t) / \{ n(r f t), 

(2.32) 



2 

3=1 &k=l j=l 



h^hi^-3a 



As usual, we shall assume that the principle of weakening spatial correlations is valid 
for the quasiequilibrium state. Then, in the thermodynamic limit, the coordinate 
distribution functions f s (r s ;t), g s (r s ;t) satisfy the boundary relations: 



(min |rj— rjfc|)->oo 



lim f s (ri,...,r s ;t) = TT nfa; t), (2.33) 

_■ — v i I i — k rv-^ 

3=1 

lim g.(ri,...,r s ;t) = 1. (2.34) 

Therefore, taking into account "slow" hydrodynamical variables (the density of the 
interaction energy in the present case) leads to a modification of boundary conditions 
for the chain of equations ( |2.29| )-( [2.31| ) for nonequilibrium distribution functions. 
In order to reproduce the usual Bogolubov boundary conditions for the weakening 
of correlations [ |22|, it is necessary to replace all g s (r s ;t) by their limiting values 
( [2.34D . Such a replacing is valid in the case of small density, however, new boundary 
conditions can be more useful for dense gasses, since they automatically take into 
account spatial correlations connected with the interaction of a pair of particles with 
the rest of the particles of the system. It is obvious that the influence of such an 
interaction increases as the density rises. 



We note that the chain of equations (|2.29| )- (|2.31| ) requires us to add equations 



for coordinate quasiequilibrium distribution functions which are functionals on the 
nonequilibrium density of the particles number n(r; t) and the inverse local tem- 
perature j3(r;t). In particular, in [ plfl it is shown that the binary quasiequilibrium 
distribution function g 2 (ri,r2;t) is connected with the pair quasiequilibrium corre- 
lation function h,2(ri,r2',t) = g2(ri,r2',t) — 1 which satisfies the Ornstein-Zernike 
equation: 

h 2 {ri, r 2 ; t) = c 2 {r x , r 2 ; t) + dr 3 c 2 {r x , r 3 ; t)h 2 (r 2 , r 3 ; t)n(r 3 ; t), (2.35) 
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where c 2 (ri,r 2 ;t) is a direct quasiequilibrium correlation function. 

Thus, the chain of equations ( [2.29| )-( j2.31| ) is distinguished from the usual Bo- 
golubov hierarchy by the existence of sources in the right-hand sides, beginning 
from the second equation, and it takes into account both one-particle and collec- 
tive hydrodynamical effects. It is important to investigate solutions to this chain 
of equations in the simplest approximations which lead to model kinetic equations 
for dense gases and simple liquids. We shall consider the most investigated binary 
collision approximation. 



2.3. Binary collision approximation 

In this section, we consider as an example of the use of hierarchy (|2.31|) , the sim- 
plest approximation of binary collisions which in the case of the Bogolubov ordinary 
conditions of correlations weakening leads to the Boltzmann equation for fi{x\\ t). In 
equation ( |2.30| ) for f 2 {x\, x 2 \ t) we omit the term with the three-particle distribution 
function, i.e. we take into account the influence of the "medium" on the evolution 
of the distinguished pair of particles only through the correlation corrections in the 
boundary condition. Then, we arrive at the equation: 



+ \L 2 + e) /2O1, x 2 ; t) = eg 2 {r u r 2 ; *)/iO&i; *)/i(»2; *)■ 



d_ 

Of 



The formal solution of (|2.36|) has the form: 



(2.36) 



u 

f 2 (x u x 2 ;t)=e j dre^+'^g^r^t + ^hix^t + ^hix^t + r). (2.37) 



Following the Abel theorem [ |9], [10], |51], |5^, [5(J, this solution can be written in the 
form: 



f 2 (xi,x 2 ;t) = lim e lL2T g 2 (r u r 2 ; t + r)/i(xi; t + r)/i(x 2 ; t + r). 



(2.38) 



Substituting expression (|2.38|) into equation (|2.29|) , one obtains a kinetic equation 
in the binary collision approximation 



d 

-Q t + ' lL i l ) ) h(x\\ t) = Ico\\{Xl, t), 



(2.39) 



where 

^coll(^l) t) 



dx 2 iL(l,2) lim e iL2T g 2 (r- 1 ,r 2 ;t + r)/ 1 (x 1 ;t + r)/ 1 (x 2 ;t + r) (2.40) 



is a collision integral. It is necessary to point out that equation fl2.35|) for the pair 
quasiequilibrium correlation function must be added to kinetic equation ( |2.29|) . 
Equation (|2.35| ) takes into account an essential part of the many-particles corre- 
lations. 
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2.3.1. Kinetic equation of the revised Enskog theory 

In papers [ [I, [| it was shown how the collision integral ( j2.40| ) in the case of a 
system of hard spheres transforms into the collision integral of the RET kinetic equa- 
tion [|34|. For subsequent manipulations, it is convenient to represent the particle 
interaction potential in the form: 

^(\r jk \) = lim n\r jk \), &{\r jk \) = { ^"Itl' (2 ' 41) 

a — >+oo I ' I jk\ ^ V) 

where a is a hard sphere diameter. Since the potential $ hs (|rjfc|) is strongly singular 
(in particular, the operator iL(l,2) is not well defined), we shall operate with the 
potential ^ a (|r*^|), when deriving the kinetic equation, and set a — > +00 in the final 
expressions. 

The limit r — > — 00 in the collision integral fl2.40| ) is mathematically formal. At 
the physical level of description, this limit assumes that |r| r , where r > is 
some characteristic time scale. Depending on the choice of To, we obtain different 
stages of the evolution of the system (kinetic, hydrodynamic) . At the kinetic stage, 
r is a characteristic interaction time for which the singular potential ( |2.41| ) is an 
arbitrarily small quantity (r — > +0). Therefore, the limit r — t — > —00 keeps 
its form even in the case r — > — 0, provided r is of the higher order of smallness 
relatively to r : 

lim J lim — I = -00. (2.42) 
In this case the collision integral (|2.40 ) [ [j], |^, Q for the hard spheres interparticle 



potential of interaction transforms into the following final expression: 

/ c h s u (x i; t) = a 2 J df 12 dv 2 6(r 12 ■ g)(r 12 ■ g) x (2.43) 

{^(r-!, n + a + f 12 ; t)/i(r 1; v[] t)/i(n + <r + r 12 , v' 2 ; t) - 
g% s {r 1 , n - a + r 12 ; t)/^, v x \ t)fi(r 1 - cx + r 12 , v 2 ; t)} , 



where 



(2.44) 



v[ = Vi + f 12 (f 12 • g), 
v' 2 = v 2 - r l2 (r l2 ■ g). 

The collision integral in the form (|2.43|) is identical to the collision integral of the 
RET theory first introduced by van Beijeren and Ernst [ 34] on the basis of a 



diagram method. As Resibois showed [ [35|, the kinetic equation of the RET 
theory satisfies an if -theorem. 

We represent the collision integral ( |2.43| ) in a more compact form by using the 
two-particle quasi-Liouvillian (evolution pseudo-operator) T(l,2) [ ^6|. Then, the 
kinetic equation takes the form: 

(Jti+Vi-^TJ = / d^ni^^ln.r^)/!^;^!^^), (2.45) 
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T(l,2) = a 2 J dr 12 6(r 12 -g)(r 12 -g)x (2.46) 
{<5(>i2 + a + r 12 )B(r 12 ) - S(r 12 - a + r 12 )} , 

B(r vt )*(v 1 ,v 2 ) = *{v' 1 ,v , 2 ), (2.47) 

where \I/ is an arbitrary function of the velocities. 

Thus, on the basis of the considered approach we have derived the kinetic equa- 
tion of the RET theory without additional phenomenological assumptions. We have 
shown that within the framework of the employed method this equation corresponds 
to the simplest pair-collision approximation without allowance for retardation in 
time. 



2.3.2. Kinetic equation for a multistep potential of interaction 

The collision integral (|2.40|) is still rather complicated to be written in an explicit 
form for an arbitrary potential of interaction. As it is now known [[30), |32j, only in two 
particular cases, namely, for hard spheres and square-well potentials, this integral is 
reduced to an analytical form. In order to extend these previous results, we consider 
the interparticle potential in the form of a multistep function 



= Hm *$(r iJk ), (2.48) 

J £()— >^ J 

where 



$Tk( r 3k) = \ e h <7,_! < r jk <<r h l = l,...,N*- (2.49) 
<7 N * < rj k ; 

and N* is the total number of walls except for a hard sphere wall. Potential of a 
hard sphere contains strong singularity (operator L(j, k) is hard to define). That is 
why we shall deal with the potential $j%('rjk) and only in final expressions we shall 
put Eq — > oo. Let us for convenience separate the systems of attractive and repulsive 
walls. Let n* be the number of repulsive walls at the distance a V i = and of the 
height Ae\ = e\ — e\ +1 > 0, where e\ = £j, % — 1, . . . , n*; and m* - the number of 
attractive walls with the corresponding parameters <j a j = cr J+n », Ae'j = e^ +1 —e^ > 0, 
where e'j = e j+n *, j — 1, . . . ,m*; e T n * +1 = e\, = 0, N* = n* + m*. Thus, the 

parameters do (hard sphere diameter), n*, <r r i, As], m*, a^j and Ae* completely 
determine the geometry of the multistep potential (see figure [I] where a specific case 
n* = 1, m* = 3 is shown and the multistep function approximates some real smooth 
potential). 

The limit r — > — oo in the collision integral (|2.40|) is a formal one in the mathe- 
matical sense. On the physical level of description this means that |r| 3> t , where 
ro > is some characteristic interval of time. For different values of To we can in- 
vestigate different stages of the evolution of the system. In the Boltzmann theory 
[ f)7| |t| ^> r c , where r c is the time of binary interactions (collision time). On the 
other hand, |r| is far less than the characteristic scale of time r m for hydrodynamic 
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8 2 



Figure 1. The multistep potential at n* = 1 and m* = 3. 

variables. The above situation is possible because of the fact that for dilute gases 
the kinetic (r ~ r c ) and hydrodynamic (r ~ r m ) stages of the evolution are far from 
one another in time, i.e. r c <C Tf <C r m , where Tf is the characteristic time of free 
motion. The pattern is qualitatively different in dense gases and liquids with a real- 
istic smooth potential of interaction. The kinetic and hydrodynamic stages appear 
to be closely connected. Besides, such quantities as the length of free motion and 
the time of interaction are not defined in a usual way, because all particles make 
influence on the dynamics of interaction for some chosen pair of particles. 

However, for special types of potentials such a classification of times remains valid 
even for high densities. For a multi-step potential the region Q of binary interactions 
consists of the following subregions [ofc— Ar , cr k + Ar ] , k = 0, 1, . . . , N*, where A 
ro — ► +0 due to the singular nature of the potential under consideration. This 
potential has the finite range max{<Tfc} = ctn* > of action. We can introduce 
the following set of specific time intervals: tq = A To/go — > +0 is the time of pair 
interactions on the walls, Ar = min{(Xfc — <Jk-i\ /go > is the time of motion 
between the two nearest neighbouring walls and t w hoie = &N*/go > is the time 
of motion of the whole systems of walls for some pair of particles, where go is an 
average relative velocity of two particles. As the potential contains the horizontal 
parts, where the force of the interparticle interaction is equal to zero, it is also 
possible to introduce time Tf as an average time of free motion of particles in the 
system. Changing the geometry of the potential and increasing the density, we can 
make this time arbitrarily small in order to match the relation 

t < r f <Ar < t who i e < r m . (2.50) 

The kinetic stage of the evolution corresponds to small times of order r ~ Tf . There- 
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fore, as far as the relation (|2.50| ) is satisfied, the formal limit r — > — oo should be 
considered as r/r — > — oo or merely as |r| ^> r . For smooth non-singular potentials 
the region of interaction has some non-zero size and To is finite. For singular poten- 
tials the case To — > +0 is possible and the limit t/tq — > — oo remains valid even for 
r — > —0 if only tq is a value of the higher infinitesimal order 



lim < lim — 

I tq-H-0 T 



-oo. 



(2.51) 



Thus, the formal limit r — > —0 can be applied to the collision integral ( |2.40| ) in the 
sense of ( 2.50 ) and (|2.51| ). As it is now well established, the limit r — > — leads to 
the kinetic equations of the RET [ |3|, |34], |38| and RDRS [ |3^j theories. 

As it was shown in [ the collision integral (|2.40 ) for potential (|2.48|) taking 
into account ( |2.5C| ), ( |2.51|) , can be presented in a compact form: 



dx 2 T 12 g 2 (ri, r 2 ; t)fi(xr, t)fi(x 2 ; t) 
in terms of the two-particle pseudo-Liouville operator 



(2.52) 



fn = f hs + Y: E ^+E E ^ 



(2.53) 



2=1 p=b,c,d 



j=l p=b,c,d 



which consists of the pseudo-Liouville operators for a hard sphere wall, for ith re- 
pulsive wall and for jth attractive wall, correspondingly: 

T hs = a 2 J df 12 (fi2 • g)6(r 12 ■ g) {5(r 12 + a+r 12 )B a {r l2 ) - 5(r 12 - a£r 12 )} , 

(2.54) 

df 12 (f 12 • g)9 p ri {r l2 ■ g) [d(r 12 + alfr 12 )B^r 12 ) - 8(r 12 - a 2 f f 12 )} , 

(2.55) 

df 12 (f 12 ■ g)9l j (r 12 • g) {5(r 12 + cr^f 12 )^.(f 12 ) - 6(r 12 - *lfr 12 )} . 

(2.56) 



T p = a 2 - 



fP. = a 2 - 



Here, 

b ri (r 1 2-g) = 0(-r u -g), 
c ri(r 12 -g) = 6{r 12 - g- al), 
6i{r l2 -g) = 6(r l2 -g)6(a\-r l2 -g), 

are the corresponding unit step functions, 



6 b aj (r 12 ■ g) = 6(r 12 ■ g) 



OlAru-g) 



0ti(ri2-g) 



1,6 
°ri 


= 


2,6 
°ri 




1,6 
aj 




2,6 

°"ai = 
a J 




l.c 
°ri 


= CT ri> 


2,c 
°ri 


= <4> 


l,c 
d J 




2,c 

°"a} = 
a J 




l,d 
°rl 


= <4> 


2,d 




X,d 




2,d 





-r 12 • g- aj), 

-^12 ■ g")0K + r l2 ■ g) 
(2.57) 



(2.58) 
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and B p are operators of the velocity displacement caused by an interaction of the 
p-type at each of the walls, namely, 

B a (r 12 )^(v u v 2 ) = Bf i {r 12 )^{v 1 ,v 2 ) = B d aj {r l2 )^{v l ,v 2 ) = *(«i,«' 2 ), 

B h ti {r 12 )y{ Vl ,v 2 ) = *« 1; < 2 ), B b aj (r 12 )*( Vl ,v 2 ) = 

B^(r 12 )^( Vl ,v 2 ) = %f 12 )*K,u 2 ) = *«i,< 2 ). 

(2.59) 



Physically, the values a r f (&' p ) (|2.58|) correspond to the distances between the 



particles just after the interaction of the p-type on the ith repulsive (jth attractive) 
wall of the potential, whereas the values o\f (erf/) correspond to the distances just 
before the collision. 

The if -theorem for the kinetic equation with a multistep potential of interaction 
( f2.52| ) was proved in [ §]. The normal solution by using the boundary conditions 



method [ fHj ^g] was found in [ [2|, . On the basis of the solution found, numerical 
calculations of transport coefficients for liquid Argon along the liquid- vapour curve 
were performed. 

2.3.3. Enskog- Landau kinetic equations for systems of charged hard 
spheres 

We consider kinetic equation ( |2.39|) in the binary collision approximation with 
collision integral ( |2.40|) , when the interaction potentials <P(\rj k \) of classical particles 



at short distances can be modelled by the hard-sphere potential <P hs (\rjk\) ( [2.41| ), 
and at large distances by a certain long-range smooth "tail" $ l {\rj k \), i.e. 

3>{\r 3k \) = ^{\r lk \) + d>\\r jk \), (2.60) 

where 

si^.n-J ' \rjk\<<r, 



We note that breaking the particle interaction potential ( |2.60|) into short- and long- 
range parts is not unique and, therefore, there arises the problem of optimal sepa- 
ration. 

These problems have often been discussed in equilibrium statistical mechanics 



[ [59|. The similar ideas about the separation of the interaction potential at the 



derivation of kinetic equations were used by Rudyak [ |60|, |6l] . 



With allowance for ( |2.3|) and ( p.41| ) , the collision integral ( p.40| ) for the interpar- 



ticle potential ( |2.60| ), ( |2.61| ) takes the form: 

Icoii(^i; t) = h{xi, t) + I 2 ( Xl ; t), (2.62) 

I x (xv,t) = - lim J dr 12 r 2 12 J dr i2 J dv 2 iL a (l,2) liir^e^ x 



o 



g 2 (n, v 2 - 1 + r)/i(xi; t + r)/i(x 2 ; t + r), (2.63) 
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h(x 2 ;t) = -J dr 12 r\ 2 j df 12 J dv 2 iL 1 (l,2)^lim o e i (^° )+Ll(1 ' 2) ) r x 

g 2 (ri, r 2 ; t + r)/i(zi; t + r)/i(x 2 ; t + r), (2.64) 

where 

Lf = m + L(2), t . (1 , 2)= ,£^ (M) (£-£) (2.65) 

and we have used the idea of time separation of the interactions: instant collisions of 
hard spheres (r — > — 0) and an extended process of interaction with the long-range 
potential (r — > — oo). In accordance with the results obtained in subsection 3.1, the 
first term in the right-hand side of (|2.62|) is identical to the collision integral of the 



hard-sphere system (the right-hand side of (|2.45| )): 



h(xi, t) = J dx 2 f (1, 2)g 2 (n, r 2 ; t)^; t)f x {x 2] t) (2.66) 

with the only difference that here £2(^1, r 2 \ t) is a quasiequilibrium binary distri- 
bution function of the particles in the system with a total interparticle interaction 
potential ( ggOfl . 

In the case when the long-range part of the interaction potential is absent 
(r) = 0), the first part I\(xx\ t) of the collision integral is identical to the collision 
integral of the RET theory; the second part I 2 (x\]t) is then identically equal to 
zero. In the second limiting case, when the density of the system is low (rarefied 
gases: n — * 0, g2i.f1, r 2 ; t) — * 1 and the hard-sphere part of the potential vanishes 
(cr — > +0), the second part I 2 {x\, t) of the collision integral is identical to the collision 



integral of the Boltzmann equation [ 30 , and at the same time Ii(x\ \ t)— >0. 

In a special case, when the long-range interaction is weak, we make an expansion 
for expKiL^ + iL^l, 2))r}, restricting ourselves to the term linear in iL^l, 2). Then, 
the second part I 2 (x\;t) ( |2.64|) reduces to the form: 

I 2 ( Xl ;t) = I 2 : °\x l ;t)+I 2 :i \x 1 ;t), (2.67) 

00 

lP(x 1 ;t) = - J dr 12 r 2 12 J df 12 J dv 2 iL l (l,2) lim e iL 2° )r x 

g 2 (ri, r 2 - 1 + r)/i(zi; t + r)h{x 2 ; t + r), (2.68) 

00 

lP( Xl ;t) = J dr l2 r\ 2 J dr 12 J dv 2 iL\l,2) lim e iL ^ T X (2.69) 
dr> e- iL( ° >r 'iL\l, 2)e lL ^ T 'g 2 (n, r 2 ; t + r)/i(zi; t + r)f x {x 2 - t + r). 



The first term I 2 °\xi, t) is a generalization of the Vlasov mean field in KMFT [ |38 
and the second I 2 (xi, t) is a generalized Landau collision integral with allowance for 
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retardation in time in the approximation of the second order on interaction. Indeed, 
if we set formally g2{r; t) = 1 (rarefied gases) and a — > 0, and completely ignore the 
spatial inhomogeneity of fi(xi,t) and the time retardation, then the second term 
I2 (xi,t) (|2.69 ) is transformed into the ordinary Landau collision integral [ |62|. In 
the case when ^ 1 (|r 12 |) is the Coulomb interaction potential, the complete collision 
integral (|2.62| ), Q2.66| ), ( |2.67|) can be called the Enskog-Landau collision integral for a 
system of charged hard spheres which, in contrast to the ordinary Landau collision 
integral, does not diverge at short distances. However, at large distances we will 
observe a divergence, as usual. To eliminate this, it is necessary to take into account 
the effects of screening [ |23J, [63, 63]]. To avoid this problem sequentially we have 
to consider the kinetic equation with taking into account the dynamical screening 
effects [ [l3j. But this way is impossible with the Enskog-Landau kinetic equation. 
The only thing we can do for further calculation is to change the upper integral limit 
to some finite value which could have a meaning of the statical screening value in 
our system (see the following subsections). To solve this problem we must consider 
dynamical screening effects. 

Using the boundary conditions method, a normal solution to the generalized 
Enskog-Landau kinetic equation was found in [ |66| . This solution coincides with 
that obtained in []3|, [7fl for a stationary case. On the basis of normal solutions, numer- 
ical calculations of such transport coefficients as viscosity and thermal conductivity 
were performed for once-ionized Argon [ [Tj and for mixtures of ionized inert gases [ 

m. 



2.4. Modified group expansions for the construction of solutions to the BBGKY 
hierarchy 

2.4.1. Modified group expansions 

Recently, a kinetic equation of the revised Enskog theory for a dense system 
of hard spheres (|2.45|) and an Enskog-Landau kinetic equation for a dense system 
of charged hard spheres with collision integrals ( |2.66| ), ( |2.67| ) have been obtained 
from the BBGKY hierarchy in the binary collisions approximation [|3j, [?]]. It should 
be noted that this approximation does not correspond to the usual two-particle 
approximation inherent in the Boltzmann theory, because an essential part of the 
many-particle correlations is implicitly taken into account by the pair quasiequilib- 
rium distribution function g2( 7 "i, r 2\ t). 

To analyze solutions to the BBGKY hierarchy ( |2.31| ) [ [D, [7] m higher approx- 
imations on interparticle correlations, it is more convenient to use the concept of 
group expansions 613] . This was applied to the BBGKY hierarchy in previous 

investigations [ ^BJ 28, |52, ^9f using the boundary conditions which correspond 
to the weakening correlations principle by Bogolubov [ p22 |. The same conception 
was envolved in papers by Zubarev and Novikov [ [52] , where a diagram method for 
obtaining solutions to the BBGKY hierarchy was developed. 

To analyze the BBGKY hierarchy fl2.31 ), we turn to the papers by Zubarev and 
Novikov [ and earlier ones by Green [ |67| and Cohen [ 65 , BIJ , and pass from the 
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nonequilibrium distribution functions f s (x s ;t) to the irreducible distribution ones 
G s (x s ;t), which can be introduced by the equalities presented in [ |52|, ^7|. In our 
case, with some modifications we obtain: 



fi(xi,t) = Gi(xi]t), 

f 2 (xi, x 2 ; t) = G 2 (x 1 , x 2 ; t) + g 2 (ri_, r 2 ; t)Gi(xi; t)G 1 (x 2 ; t), 
f 3 (xi, x 2 , x 3 ; t) = G 3 (xx,x 2 , x 3 ; t) + ^ p G 2 (xi, x 2 ; t)Gi(x 3 ; t) + 
g3(ri, r 2 , r 3 ; t)Gx{x x ] t)G x {x 2 \ t)G x {x 3 \ t), 



(2.70) 



Here, the position-dependent quasiequilibrium distribution functions g 2 (ri, r 2 ; t), 
g 3 (ri,r 2 ,r 3 ;t), g s (r s ;t) are defined in [ [IT], [T2fl . The modification of group expan- 
sions ( |2.70| ) consists in the fact that a considerable part of space time correlations 
is accumulated in the quasiequilibrium functions g s (r s ;t). If all g s (r s ;t) = 1 for 
, these group expansions coincide with those of papers [ p2l B7L BHl p9 



2,3, 



As far as each line in (|2.70| ) brings in new functions G s (r s ; t), s = 1, 2, 3, . . ., the cor- 
responding equations can be solved with respect to irreducible distribution functions 
and we can write the following: 



Gi{xi;t) = fi{xi;t), 

G 2 (x 1 ,x 2 ;t) = f 2 (xi,x 2 ;t) - g 2 (r 1} r 2 ; t)/i(xi; t)fi(x 2 ; t), 
G 3 (x 1 ,x 2 ,x 3 ;t) = f 3 {xi,x 2 ,x 3 ;t) - ^2 p f 2 (xi,x 2 ;t)f 1 (x 3 ;t) - 
h 3 (r 1 , r 2 , r 3 ; t)f\(xi, i)h(x 2 ; t)f x {x 3) t), 



(2.71) 



In (|2.70| ) and (|2.71| ), the symbol denotes the sum of all different permutations 



of coordinates for three and more particles 

h 3 (r u r 2 , r 3 ; t) = g 3 (r 1} r 2 , r 3 ; t) - g 2 (r u r 2 ; t) - g 2 (r ± , r 3 ; t) - g 2 (r 2 , r 3 ; t) 

= h' 3 ( ri ,r 2 ,r 3 ;t)-2, (2.72) 

where h' 3 (ri,r 2 ,r 3 ;t) is a three-particle quasiequilibrium correlation function. Now 
let us write the BBGKY hierarchy ( [2.3 1|) [ EE], |3L H| for the irreducible distribution 
functions G s (x s ;t), namely, the first two equations, 

J^+LL(1)J G 1 (x 1 ;t) + (2.73) 
J dx 2 iL{l,2)g 2 {r 1 ,r 2 ;t)G 1 {x 1 ;t)G 1 {x 2 ;t) + J dx 2 iL(l, 2)G 2 {x x , x 2 ; t) = 0. 



Differentiating the relation for G 2 (x\, x 2 ; t) in (|2.71|) with respect to time and using 



the second equation from the BBGKY hierarchy for the function f 2 (xi, x 2 ; t), we 
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can get for the pair irreducible distribution function G 2 (xi, x 2 ; t) an equation, which 
reads: 

d \ fd \ 

— + iL 2 + e\ G 2 (x 1 ,x 2 ;t) = - I — + iL 2 1 g 2 (r 1 ,r 2 ;t)G 1 (x 1 ;t)G 1 (x 2 ;t) - 

( dx 3 {iL(l,3) + LL(2,3)} {G 3 (x u x2,x 3 ]t) + ^G 2 {x u x 2 ;t)G 1 {x z -t) + 
J P 

g 3 (ri, r 2 , r 3 ; t)Gi(xi; t)Gi(ar 2 ; *)Gi(a; 3 ; *)}. (2.74) 

In a similar way, we can obtain other equations for the three-particle irreducible 
function G 3 (xi,x 2 , x 3 ; t) and the higher G s (x s ; t) ones. One remembers now that the 
appearance of the quasiequilibrium distribution functions g^fi, f 2 ; i), 
g 3 (r*i, r 2 , r 3 ; t), g s (r s ; t) in the hierarchy is closely connected with the fact that the 
boundary conditions for the solutions of the Liouville equation take into consider- 
ation both the nonequilibrium character of the one-particle distribution function 
and the local conservation laws, which corresponds to a consistent description of 
the kinetics and hydrodynamics of the system [ 0, 01 ■ Since in the present paper we 
analyze two first equations, fl2.73|) and Q2.74j ) only, we will not write the others. It is 
important to note that, if we put formally g s (r s ; t) = 1 for all s — 2, 3, . . . in ( [2.73D 



and ( |2.74| ), we come to the first two equations of the BBGKY hierarchy for the 
irreducible distribution functions G\{x\ \ t) and G 2 {xi, x 2 ;t), which were obtained in 
paper [ |52| by D.N.Zubarev and M.Yu.Novikov. The first term in the right-hand side 
of (|2.74j) is a peculiarity of (|2.73|) and (|2.74j ) equation system. This is a term with 



a time derivative of the pair quasiequilibrium distribution function g 2 (^i?^2;^)- As 
it was shown in [ [Lj [3| , the binary quasiequilibrium distribution function is a func- 
tional of the local values of temperature f3{r;t) and mean particle density n(r;t). 
Thus, time derivatives of g 2 (ri,r 2 \ft(t),n(t)) will conform to /3(r;t) and n(r;t). 
These quantities, in their turn, according to the self-consistency conditions [ |l|, 0, 
will be expressed via the average energy value (£'(r))* in an accompanying reference 
frame and via (n(r))', which constitute a basis of the hydrodynamical description 
of a nonequilibrium state of the system. Solving equation ( |2.74| ) for the irreducible 
quasiequilibrium two-particle distribution function G 2 (xi, x 2 \ t) in the generalized 
polarization approximation and taking into account the first equation of the chain 
fl2~73l) , lead to [0,0: 



d_ 

Of 

t 



iL(l)J G 1 (x 1 ;t) + j dx 2 iL{l,2)g 2 {r h r 2 ;t)G 1 {x 1 ;t)G 1 {x 2 ;t) = 
t 

J dx 2 J dt' e £(^ '- ^) iL(l,2)^/(M , ) x (2.75) 

—00 

(-^ + \L 2 + JSf (xi, x 2 - g 2 (n, r 2 - t')G x {x x - ti)G x (x 2 \ t'), 
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where 

U(t,t') = exp + |- jT dt" [Uv + SrixuXfX)]}, 
5? (xi, x 2 ; t) = J5f (xi; t) + JSf (x 2 ; t). (2.76) 

The operator Jtf(xi;t) can be obtained by a variation of the Vlasov collision integral 
near the nonequilibrium distribution Gi(xi,t), namely, 



dx 3 iL(l,3)G 1 (x 3 ;t)G 1 (x 1 ;t)j = (2.77) 

dx 3 LL(1, 3)Gi(x 3 ; *)5Gi(a;i; t) = Sffa; t)<SGi(zi; t). 

This is a kinetic equation for a nonequilibrium one-particle distribution function with 
the non-Markovian collision integral in the generalized polarization approximation. 
It should be noted that the presence of the Vlasov operator ^f(xi,x 2 ;t) in the 
collision integral ( |2.75| ) indicates taking into consideration collective effects. Analysis 
of the collision integral (|2.75f) in a general case is a rather complicated problem. But 
it is obvious that the collision integral in (|2.75| ) or an expression for G 2 (xi,x 2 ;t) 



with ( 2.74j) may be much simplified for every physical model of a particle system 



or for each nonequilibrium state of the collision integral in (|2.75|) . To show this, we 
shall consider two particular cases: a hard spheres model and Coulomb plasma. 

2.4.2. Hard spheres model in the polarization approximation 

In this subsection we shall investigate kinetic processes for a hard spheres model 
in the approximations which are higher than a binary collisions one. We take into 
account the character of the model parameters and the results of the previous section 
of this article and papers [ [H| |49| . This investigation is convenient to carry out on the 
basis of the equation chain fl2.73|) , ( 2.74Q at the formal substitution of a potential 



part of the Liouville operator iL(l, 2) by the Enskog collision operator T(l, 2) [ 
26], ITOf. In this case, equations (|2.73|) and (|2.74 ) have the form: 



^ + iL(l)^G 1 (x 1 ;t)+ (2.78) 
J dx 2 f{l,2)g 2 {r 1 ,r 2 ;t)G 1 (x 1 ;t)G 1 (x 2 ;t) + J dx 2 f (1, 2)G 2 ( Xl , x 2 ; t) = 0, 

^ + \L° 2 + f (1, 2) + e^j G 2 { Xl , x 2 - 1) = (2.79) 

- I d.x 3 {f(l,3)+f(2,3)}{G 3 (x 1 ,a; 2 ,X3;t) + ^G 2 (x 1 ,a; 2 ;t)G 1 (x3;t) 
J P 

+ g3(n, r 2 , r 3 ; t)G?i(xi; t)Gi(x 2 ; *)C?i(x 3 ; t) \. 



+ iL° 2 + f (1, 2) ) g 2 (n, r 2 ; i)^^; t)G x (x* t) 
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Further, we will consider the same approximations concerning equation Q2.79Q , in 
which G 3 (xi,x 2 , x 3 ; t) and hs(ri, r 2 , r 3 ; t) are neglected. Then, if we introduce simi- 
larly to (|2.77|) the Boltzmann-Enskog collision operator C{x\\ t) (££{x\\ t)—*C{x\ \ £)), 
equation ( |2.79| ) could be rewritten in the next form: 



5 J dx 3 T(l, 3)Gi(xi; t)Gi(z 3 ; t) = C{x x - t)6G 1 (x 1 ; t), (2.80) 
(jt + iL° 2 + f (1, 2) + C(x u x 2 ; t) + e^j G 2 (x 1 , x 2 ; t) = (2.81) 
- (§1 + iL °2 + T(l, 2) + C(xi, x 2 - 1)\ g 2 (n, r 2 ; t)G x { Xl \ t)G x {x 2 - 1), 
Hence it appears that the formal solution to G 2 (xi,x 2 ;t) reads: 



G 2 {x 1 , x 2 ; t) = - J dt' e^'-^U hs (t, x (2.82) 

— oo 

l ^ + \L\ + T(l, 2) + C(x 1; x 2 ; t') | g 2 (r 1; r 2 ; Od^i; O^iO^ 0, 
where Uhs(t,t') is an evolution operator for the system of hard spheres: 



U hs (t,t') = exp_ 



dt" 



iL^ + ftl^ + COr!,^") 



(2.83) 



C(x 1 ,x 2 ; t) = C(xi\ t) + C(x 2 ] t). 



Now let us put (|2.82j ) into the first equation ( |2.78|) . Then, the resulting equation 
takes the form 

^+\L{l) S jG 1 {x 1 ;t) = J dx 2 f(l,2)g 2 {r 1 ,r 2 ;t)G 1 (x 1 ;t)G 1 (x 2 ;t)- (2.84) 
o 

dx 2 T(l,2) J dt' e^'-VUUt,t'){^ + iL° + f (1,2) + C( Xl ,x 2 ;t')^ x 

— oo 

g 2 (r 1 ,r 2 ;t')G 1 (x 1 ;t')G 1 (x 2 ;t'). 

This equation can be called a generalized kinetic equation for the nonequilibrium 
one-particle distribution function of hard spheres with a non-Markovian collision 
integral in the generalized polarization approximation. The first term in the right- 
hand side of this equation is the collision integral from the revised Enskog theory 
[ U, [7|, |34], |54j| . Neglecting time retardation effects and assuming that the operator 
C(xi,x 2 ;t) does not depend on time when 

G 1 (x 1 ; t) = /o (p) = n (^) ^exp {~^} 
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is a local equilibrium Maxwell distribution function, the next term can be rewritten 
in a simplified form: 

t t 
I R (x 1 ;t) = - J dt' e^'-^Roix^t, t')G x {x X ] t') - J dt' e^^R^-t^G^t'), 

— oo — oo 

(2.85) 

where 

R (x x ;t,t') = J ds 2 f(l,2)exp|(t'-t) [iL° + T(l, 2) + C(x x , x 2 )] } x 

[iL° + C(x x ,x 2 )] g 2 (r 1 , r 2 ; t')G x (x 2 ; t'), (2.86) 

R x {x x ;t,t') = y"dx 2 f(l,2)exp|(t'-t) iL° + f(l, 2) + Cfo, ar 2 ) }x 

fCl^JgaCrLraiOGiCxajO, (2-87) 

t, t') is a generalized ring operator. The kinetic equation Q2.84p together with 
(p.86|) and (|2.87 ) is the generalization of the kinetic equation for a system of hard 



spheres which was obtained by Bogolubov in [ |49|, |70[. It coincides with the case, 
when the quasiequilibrium pair distribution function of hard spheres is set formally 
to be a unity. 

2.4.3. Coulomb plasma in the polarization approximation 

Here we shall study an electron gas, which is contained in a homogeneous pos- 
itively charged equilibrating background. This background can be created, for ex- 
ample, by hard motionless ions. Then, electrons interact according to the Coulomb 
law: 



<P(\r 



12 



\r x -r 2 \ \r 12 \ 

the Fourier transform of which exists in the form of a real function <?(|fc|): 



e 

r X2 



dk 



■*(|fc|)e 



ifcri 



<2>(|fc|) = <P(k) 



-2 ' 



{21 



here k is a wavevector, e is an electron charge. Let us consider equation chain (|2.73 ), 
( P-74|) , when G 3 (x x , x 2 , x 3 ; t) = 0, g 3 (ri, r 2 , r 3 ; t) — in the homogeneous case, when 
G x (x x ;t) = Gi(pi,t) and pair distribution functions depend on |ri 2 |. Following the 
Bogolubov method [ p2| , we shall assume that the one-particle distribution func- 
tion Gi(pi,t) is calculated in the "zeroth" order on the interaction constant q, pair 
distribution functions G 2 (r X2 ,p 1 ,p 2 ; t) and g 2 (r X2 ;t) in the first order on q, and 
G 3 (xi,x 2 , x 3 ; t), g 3 {r x , r 2 , r 3 ; t) ~ q 2 , where q = ^0, r d = ^0/Ane 2 n is the Debye 
radius, n = N/V, = k^T, ks is the Boltzmann constant, T is thermodynamic 
temperature. Therefore, to obtain an equation for G 2 (r X2 ,p 1 ,p 2 ; t) in the first ap- 
proximation on the interaction constant q without time retardment effects, it is 



26 



A consistent description of kinetics and hydrodynamics 



necessary to retain all integral terms, but omit all the others. In this case, using the 
Fourier transform with respect to spatial coordinates for a homogeneous Coulomb 
electron gas, the set of equations (|2.73Q , ( j2.74| ) yields: 



g- t Gi( Pl ;t) = -— J dkdp 2 i^(ife|)g 2 (fc;t)G 1 (p 1 ;t)G 1 (p 2 ;t) 

d r 

dkdp 2 i<P(\k\)G 2 (k, Pl ,p 2 ;t) : 



dpi 



or 



Q-Gi{p^t) = —G 1 (p 1 ;t) J dk k$(\k\)Zmg 2 (k;t) 

+ [ dkk$(\k\)%mG 2 (k lPl ;t) (2.89) 

and an equation for G 2 {k,p 1) p 2) t): 

^- + ik^+e^G 2 (k, Pl ,p 2 ;t) = (2.90) 

4+ 



ik$(\k\) l-J^-Gi( Pl -,t) J dp 3 G 2 (k,p 2 ,p 3 ;t) - -^-G 1 (p 2 ;t) J dp 3 G 2 (k, Pl ,p 3 ; 
ifep(|fc|) | — G ^ Pl]t ) g2 {-k- t)G x (p 2 - 1) - —Gife; % 2 (fe; t)Gi(p l5 1) 



e —> +0, and G 2 (k, Pl ;t) = j dp 2 G 2 (k,p 1 ,p 2 ;t); Qmg 2 (k;t), QmG 2 (xi,x 2 ;t) are 
imaginary parts of the corresponding distribution functions. The following properties 
should be noted: 

G 2 (-k, Pl ,p 2 ;t) = GK^p^p^t), 
g 2 (-k;t) = g* 2 (k;t), 



where * denotes a complex conjugation. The solution to ( |2.90|) , neglecting time 
retardment effects, reads: 

G 2 (k lPl ,p 2 ;t)= (2.91) 

k<P(\k\) (8 8 
' M 1 -G 1 (p 1 ;t)G 2 (-k,p 2 ;t)-—G 1 (p 2 ;t)G 2 (k,p 1 ;t) 



k.Pn-io\d Pl ^ J ZK d P2 



m 



k<P(\k\) (8 8 

11 1 -G 1 (p 1 -t)g 2 (-k ] t)G 1 (p 2] t)-—G 1 (p 2] t)g 2 (k ] t)G 1 ( Pl -^ 



fc.^-iol^i ^ '~ V 8p 2 



m 



It should also be noted that equation Q2.89| ) contains an imaginary part of the 



irreducible pair nonequilibrium distribution function, to be integrated with respect 
to momentum of the second particle. Now one integrates equation ( |2.91| ) over all the 
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values of momentum p 2 and defines in such a way some function G 2 {k,p 1 ; t): 



k<P(\k\) d „ , 



k .Pll- i0 d P2 

m 



a 



G 2 (k, Pl ;t)= (2.92) 

^-(h(pi;t) [ dp 2 k pl{ k[} G 2 (-k,p 2 ;t) + 

°Px J k ■ ^ - iO 

m 

j dp 2 k p^ 2 kl \ |^-G 1 (p 1 ;% 2 (-fe;i)G 1 (p 2 ;t) - ^^(p,; t)g 2 (k; tJGifo; i)| . 
m 

Further, we should exclude from (|2.92| ) the term with G 2 (—k,p 2 ;t). To do this, we 



follow Lenard [ 47, [|!| and integrate equation ( 2.92 ) over the momentum component 
Pi±, which is perpendicular to wavevector k. The resulting expression then reads: 

8 f k<P(\k\) 

[l + $(\k\)X{k,pi;t)]G 2 (k, Pl ;t) = —G 1 (p 1 ;t) / dp 2 » G 2 (-k,p 2 ;t) + 

5 Pi ■/ • — - iO 

m 

k&(\k\) (8 8 
dp 2 prl 1 —G 1 ( Pl ;t)g 2 (-k;t)G 1 (p 2 ;t) - — G 1 (p 2 ;% 2 (fe;t)G 1 (p 1 ;t) ,. 



77? 



(2.93) 



Here the following conventional designations have been introduced: 



X(k,pi,t) = / dp 2 



k 







m 

Gi(pr,t) = I dp 1± Giip^t), 
G 2 (k, Pl ;t) = J dp 1± G 2 (k,p 1 \t) 



Gi(p 2 ;t), pi 



P'2 



Pi ■ k 

k 

p 2 ■ k 



(2.94) 



k 

k = \k\. 



Now we multiply the both equations, ( 2.92|) and ( |2.93|) , by -^Gi(pi,t) and by 



YjGi(pi,t), respectively, and subtract them: 
(l + ^(\k\)X(k, Pl ;t) 
<P(\k\)X(k, Pl ;t)g 2 (k;t) 



d d 
G 2 (fe,p 1 ;t)— Gi(pi;t) - G 2 (k,p 1 ;t)^—G 1 (p 1 ;t) 
opi opi 



G 1 (p 1 ; t) —G 1 {p 1 ;t)-G 1 [p x ; t) —G, ( Pl ; t) 



. (2.95) 



If we extract the imaginary part of this equation, one can find the unknown quantity 
$5mG 2 (fc,p 1 ; t), provided ^smG 2 (k, Pl ;t) = [|TT|, 



_d_ 

dpi 



■Gi(pi;t) 9fmG f 2 (fc,p 1 ;t) 



gjfcQgm [X(fc,pi;t)g 2 (M)] 
|l+^(|fc|)X(fc, Pl ;t)| 2 



x 



G!(p i; t)^-Gi(p i; *) - Gi(p i; t) «-Gi(pi; t) 
o P i opi 



(2.96) 



28 



A consistent description of kinetics and hydrodynamics 



Since ^smX(k,pi;t)=-ir-^-Gi(pi;t) \ |64l , putting an expression for 5mG2(fc,p 1 ;t) 
into equation (|2.89|) gives the generalized Bogolubov-Lenard-Balescu kinetic equa- 
tion for an electron gas in an equilibrating background: 



(2.97) 



d 

9pi 



dp 2 Q(pi,p 2 ;*) 



d 



_dp 1 dp 2 



where Q(Px,P2',t) is a second rank tensor 



Q{Pi,P 2 ,t) = -7T / dfc 



\<P(\k\)\ 2 k ■ k 



l + 0(\k\)X(k, Pl ;t)f 



■%mg 2 (k-t)5(k-( Pl -p 2 )), (2.98) 



which coincides with Q(p 1 ,p 2 ) [ at ^sxng 2 (k; t) — 1. In this case, the kinetic equa- 
tion (|2.97| ) transforms into the well-known Lenard-Balescu equation [ [IT], ^9|. |64 



Evidently, the generalized Bogolubov-Lenard-Balescu kinetic equation fl2.97Q claims 
the description of a dense electron gas, since in both the generalized mean field and 
the generalized Bogolubov-Lenard-Balescu collision integrals, many-particle corre- 
lations are treated by the imaginary part of g 2 (k;t). Nevertheless, the problem of 
divergence in the collision integral of equation ( |2.97| ) at small distances (k — > oo) 
still remains. There are papers where the divergence of collision integrals is avoided 
with the help of a special choice of the differential cross section (quantum systems 
[ ED])) or y i a a combination of simpler collision integrals (classical systems [ |24j1). 
These generalizations for collision integrals are attractive by their simplicity and 
helpful for ideal plasma. But, contrary to the obtained by us Bogolubov-Lenard- 
Balescu kinetic equation, they do not work for nonideal plasma. In accordance with 
the proposed structure of the collision integral 



'total 



'Boltzmann 



' Landau 



Lenard— Balescu 



the influence of particles interaction on plasma energy will be defined by the corre- 
lation function g 2 (r). Its asymptotic is 



lim 

r— >oo 



exp 



e a e b 



Boltzmann 



1 + 



e a e b 



e a e b c - r /r D 



rk^T rkftT 

Landau Lenard— Balescu. 



where r D denotes the Debye radius. In other combinations one arrives at false ex- 
pressions for thermodynamic functions [ 24] . Dynamical screening, which appears 
in the generalized Bogolubov-Lenard-Balescu collision integral obtained by us, is 
free of these discrepancies. Generally speaking, the problem of divergency could be 
solved within the framework of a charged hard spheres model, combining the results 
of this section and the preciding one. But this step is an intricate and complicated 
problem and needs a separate consideration. 
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Evidently, an investigation of the obtained kinetic equation is important in view 
of its solutions and studying transport coefficients and time correlation functions for 
model systems. 

In view of the dense systems study, where the consideration of spatial interparti- 
cle correlations is important, the BBGKY hierarchy ( 2.73Q , ( |2 . 74 ) with the modified 
boundary conditions and group expansions has quite a good perspective. The ki- 
netic equation (&M )- (&87\) is a generalization of the Bogolubov one [ f|9|, [70j for a 



system of hard spheres. M.Ernst and J.Dorfman [ [72]] investigated collective modes 
in an inhomogeneous gas and showed that the solution of a dispersion equation for 
hydrodynamic modes leads to a nonanalytic frequency dependence on a wavevector. 
This is connected with the fact that the ring operator for inhomogeneous systems at 
small wavenumbers has a term proportional to y/k. Similar investigations of collec- 
tive modes and time correlation functions in the hydrodynamic region were carried 
out by Bogolubov [ |70|. Nevertheless, it is necessary to carry out analogous inves- 
tigations of hydrodynamic collective modes and time correlation functions on the 
basis of kinetic equation ( |2.84| ), taking into account ( |2.85| )-( |2.87| ), where some part 
of space correlations is considered in the pair quasiequilibrium function £2(^1? r 2',t). 
Obviously, these results may appear to be good for very dense gases, which could 
be described by a hard spheres model. An important factor is that in the kinetic 
equation (pMD - flTSTD , as well as in the generalized Bogolubov-Lenard-Balescu one, 
collective effects are taken into account both via the Vlasov mean field and the bi- 
nary quasiequilibrium correlation function which is a functional of nonequilibrium 
values of temperature and a chemical potential. 

Transferring the obtained results to quantum systems is not obvious. Such a 
procedure is rather complicated and needs additional investigations. Nevertheless, 



some steps in this way have been already done by Morozov and Ropke [ 0, [TO, [16 



3. Generalized transport equations and time correlation 
functions 



3.1. Overview 

One of the main problems of the nonequilibrium statistical mechanics of liquids 
is an investigation of collective excitations, time correlation functions and transport 
coefficients because, using these quantities, we can compare the corresponding theo- 
retical results with the data on light and neutron scattering as well as with the results 
obtained by the method of molecular dynamics (MD) [ [73], [74], [75], [76], |77|, [78], |79 



At present, one can distinguish three regions in these investigations. 

The first one is the usual hydrodynamics in a linear approximation when the 
time and spatial evolution of small excitations in liquids can be described in terms 
of hydrodynamic modes (in terms of eigenvalues and eigenvectors of the linear hy- 
drodynamics equations) 



31, 77, 80, 31, 82 



the heat mode z^(k) 
viscosity modes z Ul 2 - 



26], [31]) HZl HP]) §2|. Such hydrodynamic modes are 
-D^k 2 , two sound modes z±(k) = ±ick — Tk 2 and two 
— vk 2 (see, for example, f |8TH). In our notations k denotes 
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a wavenumber which describes spatial dispersion of excitations, = X/nc p is 
a coefficient of thermal diffusion, v = rj/mn is the kinematical viscosity, T = 
\y + ~K/mn + \{c v jcy — l)D T is the constant of sound decrement. A, k and 77 
are the thermal conductivity, bulk and shear viscosity coefficients defined by the 
Green-Kubo formulas, n is the average density, m is the mass of a separate par- 
ticle, c p and cy are specific heats at constant pressure and volume, c is the sound 
velocity. In this region, where is much larger than atomic sizes of liquid a, 

and ujt <C 1, where uj is frequency and r is the characteristic correlation time, the 
dynamical structure factor S(k\ uj) is well described by the Landau-Plachek formula 
[ |26|, [H], |84|] in terms of Brillouin lines which corresponds to heat and sound modes. 
With increasing wavevector values k, the precision of the hydrodynamic description 
decreases, since correlations at short times and small spatial distances which are 
inherent in neutron scattering in liquids are not described by linear hydrodynamic 
equations. 

Extensions of the usual hydrodynamics were performed on the basis of modern 
methods of nonequilibrium statistical mechanics by Zwanzig and Mori [ |85|, |86|, |87 



with the use of the method of projection operators, by Sergeev and Tishchenko 



B9[ using the NSO method [ 50, 51] and by Tserkovnikov with the use of the 
method of Green functions [ |90, [H], |92|. Transport equations were obtained for 



the average values of the densities of mass, momentum and energy. These equations 
generalize the hydrodynamic ones. In these equations, thermodynamic quantities c p , 
cy depend on wavevector k, whereas the transport coefficients A, k and rj depend on 
k, as well as on frequency uj. Moreover, new generalized transport coefficients appear 
in the generalized hydrodynamic equations. These coefficients describe dynamical 
correlations between the thermal and viscous motions which vanish in the limit 
k — > and uj — > 0. This is the second region, namely, the region of the generalized 

Papers by Gotze, Liicke, 



hydrodynamics [ ||, [77], |7|, gjjj [89], ||, |91], ||, |9S 

Bosse and others [ 95, 96, 97, 98, 9£, 100 , 101 ] are of special interest. In these 



papers, an investigation of spectra for fluctuations of densities for the number of 
particles, their longitudinal and transverse fluxes for small as well as for intermediate 



values of k and uj is performed. For Argon [ [73], [76j and Rubidium [ [73], |75|, |76 



a good coincidence with the experimental data was obtained. It is necessary to 

[lOl HOR [1061 



point out papers by Yulmetyev and Shurygin [ |102| , |103 



1Q3| , [104 , [105| and Argon 



where the 

m 



dynamical structure factor S(k; uj) for Rubidium 
|106l was investigated by the projection operator method taking into account non 
Markovian effects. The results obtained in these papers are valid in the low-frequency 
approximation. 

The third region of investigations on collective modes and time correlation func- 
tions in dense gases and liquids is connected with the kinetic theory [ [77], [LOT], |108| . 
on the basis of the method of projection operators by Mori and its general- 
In this approach, the nonequilibrium one-particle distribution 



109 



9 B7 



izations 

function in phase space of coordinates and momenta is a variable for an abbreviated 
description of a nonequilibrium state of the system. At the same time, collective 
effects arise explicitly in memory functions. Approximate calculations of memory 
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functions (expansions on density, weak interaction) in the hydrodynamic limit were 
performed in papers by Mazenko [ |109| , |110| , |111| , |1 12| , |113| , |1 14|| , Forster and Martin 
j |108|| , Forster [ I15|| and others [ 100 , 116 , 117 1. John and Forster unified paper 



[ |118|| and proposed a formalism similar to that of the generalized hydrodynam- 
ics [ [103, |U|, [109], |11|, pl^l . This formalism is distinguished by the fact that the 
nonequilibrium one-particle distribution function in phase space of coordinates and 
momenta, together with the density of total energy are included into a set of variables 
of an abbreviated description. The results of this theory for the dynamical structure 
factor S(k; to) of Argon agree well with the experimental data and MD calculations 
in the regions of small and intermediate values of k and to. The same formalism was 
used later in papers by Sjodin and Sjolander [ |119|| in which the dynamical structure 
factor S(k; uj) for Rubidium was investigated. The authors of this paper touched 
a number of interesting questions about the influence of one-particle motion and 
temperature fluctuations on S(k; uj) in the regions of intermediate and great values 
of wavenumber k, as well as about the existence of short-length collective modes [ 



75]. Such questions and a number of others arise in discussing the results of neutron 



scattering with I ~ a, (|fe| = 2tt/£), where the kinetic theory of liquids is developed 



insufficiently From the experimental results of neutron scattering in Rubidium [ |75| 
it follows that there are sharp side peaks for S(k; uj) at points uj = ±ck, where c is 
an adiabatic velocity of sound. Such side peaks could be identified as Brillouin ones. 
However, the Landau-Plachek formula is not valid in this region, where I ~ a. 



During the last period 
modes, beginning from the papers by de Schepper and Cohen 
sively discussed 



the question of the existence of short-length collective 

121| , is inten- 



120 



The proposed in [ |120| , |121|] approach is based on the kinetic equation of the revised 
Enskog theory (RET) [ p^j , which was derived by us sequentially in subsection 
2.3.1. The linearized kinetic equation of RET for hard spheres, which was used in 
T201 



121 



was obtained in papers by Mazenko and others as well [ lT2", 113 , 117 
Although the kinetic theory considered by de Schepper and Cohen for short-length 
collective modes is related to a model system of hard spheres, it has explained many 
questions which are significant for understanding collective modes in the region of 
large k. First of all, it proposed a concept of extending hydrodynamic modes into 
the region of large wavevectors [ |121| , |122j| . Solving the problem on eigenvalues for 
the generalized Enskog operator (taking into account a mean field) for hard spheres 
has shown [ |122| , |132|1 the existence of five generalized (extended) hydrodynamic and 



kinetic modes. The generalized hydrodynamic modes are a generalization of hydro- 
dynamic modes to the region of large k and tend to zero at k — > 0. The kinetic modes 
at k — > take positive nonzero values. The dynamical structure factor S(k; uj) is 
presented in this approach in the form of a sum of spectral terms which correspond 
to the Lorentzian form of lines. Secondly, to confirm the theoretical predictions of 
papers [ [120L [mj |l22l |123], gH , [126| , |T27| , [128| , [1291 , [130| , [131] , [132|, [133|, [Eyj, H , 
experiments on neutron scattering in liquid Argon [ |123j| , Neon [ |125|| , as well as MD 
calculations for a system of hard spheres [ 126] and for a liquid with the Lennard- 
Jones-like interparticle potential of interaction [ |128|| were made. In [ |123|| short- 
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length sound modes and their decrement in the region 
were investigated. The obtained results 



a for a liquid Argon 



1 3 1|| for a nonanalytical dependence of the 



sound dispersion on k were compared with the experimental data and the results of 
the theory of coupled modes [ |136|| . Eigenvalues for the heat mode z^(k), obtained 
on the basis of the generalized Enskog equation, were investigated in [ |130|| . The 
calculations performed agree quantitatively with the MD data for hard spheres [ 



125 1 and with the data on neutron scattering for Argon [ 123, 124, 137] and Crypton 



138 



Interesting results of the investigation of the dynamical structure factor S(k; uo) 
and time correlation functions of flow particle and enthalpy densities, as well as 
collective modes for a liquid with the Lennard- Jones potential of interaction were 
obtained in [ 139|| . There a system of equations for time correlation functions of den- 



sities for the number of particles, momentum and enthalpy, the generalized stress 
tensor and the flow of energy obtained by the Mori projection operator method 
was used. The basis of this system are equations of a generalized hydrodynamic 
description obtained for the first time in [ |93| for simple liquid and ionic systems 
with the use of the NSO method 



140 



Solving the system of equations for time 
correlation functions in the Markovian approximation on eigenvalues in the hydro- 
dynamic limit k — > 0, uo — > showed [ 139 1 that apart from the eigenvalues which 



correspond to pure hydrodynamic modes, in particular, to heat z^{k) and two sound 
modes z±(k), there exist two eigenvalues that correspond to kinetic modes. They 
differ from zero in the limit k — > 0. The analogous results were obtained in our 
paper [ |141|| . Further, such an approach was developed for a Lennard- Jones fluid in 

TO EM IH 



a number of works [ |4|, [143], £3|, |4|, [146], [14 



151] , [152| 153, [154] 



The concept of generalized collective modes was developed [ [1421 |143|] on the ba- 
sis of equations of extended hydrodynamics and the relevant equilibrium static 
quantities which normally are obtained within the method of molecular dynam- 
ics (MD). This approach gives the possibility to calculate generalized-mode spectra 
of a Lennard- Jones fluid using nine- and four-mode descriptions for longitudinal [ 
142], |147|| and transverse [ |142| , |146j| fluctuations, respectively. The investigation of 



generalized transport coefficients dependent on wavevector and frequency was car- 
ried out for the first time in [ |148|| . Further, this approach was developed for mixtures 



water 



154] , |155| , |156|| and polar liquids, in particular, for Stockmayer models and TIP4P 
I57|, |I5B|, [ToD], [IM], [T5g] ggg,[TBJ, [Teg. 
Problems of the description of kinetic and hydrodynamic fluctuations, taking 
into account correlations in a simple fluid, were discussed in papers by Peletminskii, 
Sokolovskii and Slusarenko [ |166| , p.67|| . They used a functional hypothesis and the 
method of an abbreviated description [ |168|| . 

On the basis of general equations of fluctuational hydrodynamics obtained in 
[ 166| , we carried out investigations of liquid hydrodynamics near the equilibrium 
without taking into account far correlations and spatial dispersion of transport co- 
efficients [[1691 TO- 



It is necessary to point out a series of papers by Balabanyan [ |171|, 172 , 173 , 174 , 



175|| devoted to the investigation of classical time correlation and Poisson Green 
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functions and spectra of collective modes on the basis of the Boltzmann kinetic 
equation by the method of moments for a system of hard spheres and Maxwellian 
molecules. 

The problem of the construction of interpolation formulas for density- density 
and flow-flow correlation functions which are valid for a wide range of frequencies 
and wavevectors was considered in papers by Tserkovnikov [ fll], f^, |176| , |177|1 . Here a 
method of Green functions was applied to the molecular hydrodynamics of quantum 
Bose-systems. These questions were presented in detail in [ 176 , 177 1 for a weakly 



nonideal Bose-gas. The obtained interpolation formulas for the Green functions of 
transverse components of flow density [ |176|| , as well as of fluctuations of the number 
of particles and energy [ |177|| , appear to be valid in the hydrodynamic region and in 
the region of great frequencies and small wavelengths. These results are obviously 
valid in the classical well. 

From the presented survey, one can draw the following conclusions: firstly, for a 
more detailed study of collective modes, generalized transport coefficients and time 
correlation functions in liquids and dense gases, it is necessary to have a theory in 
which kinetic and hydrodynamic processes are considered simultaneously; secondly, 
at the present time, a consistent theory for the description of light and neutron scat- 
tering in the whole range of wavevector k and frequency u has not been formulated 
yet. 

In our papers [ U EE & 01 we P ro P ose d an approach in which the kinetics and 
hydrodynamics of transport processes in dense gases and liquids are considered as 
coupled. As a result, we obtained a system of coupled generalized equations for the 
nonequilibrium one-particle distribution function and the mean density of the total 
energy. It can be applied to the description of nonequilibrium states of the system 
of particles which are as close as they are far from the equilibrium state. In the next 
subsections of this section we shall present this approach in comparison with other 
theories and investigate time correlation functions, generalized transport coefficients 
and collective modes in comparison with the molecular hydrodynamics [ [T7|] and its 
extended versions 



139, 141, 142, 143, 144, MB, 146, 147 



149, 150, 151, 152] 



3.2. The generalized kinetic equation for the nonequilibrium one-particle dis- 
tribution function with taking into account transport equations for the 
mean density of energy 

According to the formulation of a modified boundary condition for the Liouville 
equation ( |2.4j ) with choosing parameters of an abbreviated description (ni(x)Y = 
fi(x;t) and (S(r)Y (the definition of the quasiequilibrium distribution function is 
presented by equation ( |2.12| )), the dependence of the the quasiequilibrium distribu- 
tion function on time is defined by variation in time of the average values (ni(x)}* 
and {E{r)Y: 

Q {x N ] t)=g( K ... ) {{h 1 {x))\{£{r)Y}..). (3.1) 

Time evolution of the parameters of an abbreviated description f\(x;t), (S(r)Y is 
described by transport equations. For obtaining these equations it is more convenient 
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to use the method of projection operators, which was widely used by Robertson [ 
178| , [L79|| and modified by Kawasaki and Gunton [ |180|| . We shall reformulate the 



NSO method in order to take into account projection operators. 
Let us write the Liouville equation with source ( |2.4| ) in the form: 



(J^ + iL N + e^j Ag (x N ; t) = - + iL N ^j g q (x N ; t) , 



(3.2) 



introducing Ag[x N ;t) = g(x N ;t) — g q (x N ;t). We shall use the quasiequilibrium 
distribution function g q (x N ;i) defined as follows: 

£ q (x N ;t) = exp - J dr (3(r;t)£(r) - J dx b(x; i)rii(x)j , (3.3) 



which is found from ( [2.12| ), taking into account (|2.13|) , by redefinition of the param- 



eter conjugated to (ni(x)Y: 

2 

b{x-t) = a{x-t)-(3{r-t)?-. (3.4) 

2m 

Parameters /3(r;t) and b(x;t) in ( |3.3D are defined from the corresponding self- 
consistency conditions: 

<£(r)>* = <£(r)>«, (3 5) 

(ni(x)Y = (ni(a;))*. 

Taking into account the structure of the quasiequilibrium distribution function 

m 



]3|), the time derivative 4 of this function in the right-hand side of equation ( [3^2] 



can be presented as 

J^q {x N ; t) = -g> q (t)iL N Q (x N ; t) , (3.6) 
where & q (t) is the Kawasaki- Gunton projection operator, defined as follows: 

»4t)d = (3.7) 

* {x "-' t] -J dr im?^' - 1 dx / dr ~ »' + 

This operator acts only on distribution functions and has the properties: 

& q (t)g(x N ;t') = g q (x N ;t), 
&> q (t)g q (x N ;t>) = g q {x N ;t), 



q 



(*W) = & q (t). 
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Taking into account ( |3.6|) , the Liouville equation ( |2.4|) transforms into the form 
(J- t - (l-P^iLN + e) Ag(x N ;t) = -(l - ^ q (t)yL N0q (x N ;t) , (3.8) 

the formal solution to which is 

t 

Ag (x N ; t) = - J dt' e e( ''"* } T(t, t') (l - ^ q (t')) iL N p q (x N ; t') , (3.9) 

— oo 

and one finds an expression for the nonequilibrium distribution function 

t 

Q {x N ;t) = g q (x N ;t) - J dt' e< t '-*)T(t,t')(l - & > q (t'))iL N g q (x N ;f) , (3.10) 

— oo 

where t 

T(M') = ex P+ {- jf dt" (l - ^ q (0)i£;v} (3.11) 

is a generalized evolution operator with taking account the Kawasaki- Gunton projec- 
tion operator, exp + is an ordering exponent. Solution ( |3.10|) is exact. It corresponds 



to the idea of an abbreviated description of a nonequilibrium state of the system 
( [3.1D . Acting by the operators (1 — & q (t')) and iL N on g q (x N ;t') in the right-hand 
side of expression ( |3.10|) , we write an expression for the nonequilibrium distribution 



function in the explicit form: 

g{x N ;t)=g q (x N ;t) + [ dr I dt' e e ^T{t, t')I £ (r, f)/3(r; t') 6q (x N ; t') 



-oo 

t 



dx / dt' e^'-^Tit, t')I n {x, t')b(x; t')g q (x N ; t') , (3.12) 



where 

I £ (r,t') = (l-&(t'))£(r), (3.13) 

I n (x,t') = (l-^(O)Ai(a:) (3.14) 

are generalized flows, 8{r) = iL]y£(r), hi{x) = iLjvni(x). &{t') is a generalized 
Mori projection operator which acts on dynamical variables A(r) and has the fol- 
lowing structure: 

0?(t)A{r') = <i(r')>* + (3.15) 



d(£(r)Y 



(£{r) - (S(r)Y) + jdx ±±jk ( M x) - («!(*)>*) . 
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The projection operator &{t) has the following properties: 

&>(t)£(r) = £{r), &>(t)(l - &(t)) = 0, 

@>{t)hi{x) = n x {x), &>{t)@>{t!) = ^(t). 



(3.16) 



With the help of the solution to the Liouville equation (p. 12 ), one obtains a system of 
coupled equations for the nonequilibrium one-particle distribution function fi(x;t) 
and the average density of energy (£(r))*. To do it explicitly, it is necessary to 
calculate 



Of 



{S(r)Y = {£{r))\ 



To this end let us use the equalities 



l-^(t)Ui(x)) = - (hi(x)) 



l-&(t))i(r)) = (i(r))* - (£(r))l 



and perform averaging in the right-hand side with the nonequilibrium distribution 
function ( |3.12| ). Then we obtain, calculating (ni(x))q, (£(r))q [ |[: 



d_ 

at 



fi(x;t) + ^Mx;t) = J dx' ^(\ r -r'\)^g 2 (r,r , ;t)f 1 (x;t)f 1 (x'- 1 t) + 

t 

dx' [ dt' e £(t '- t) cp nn (x,x';t,t')b{x';t') + 



— oo 

f 



dr' / d«'e ! <''-'W(zy ; «, (3.17) 



dr' dp dp' -j-<P(\r - r'|)g 2 (r,r / ;t)/ 1 (x;0A(^ , ;t) (3.18) 
m or 



dr' dp dp' — — 

m or 



, p d \ p 2 I 



— + -<P(\r-r'\) 
2m 2 Vl u 



g 2 {r,r'-t)h{x-t)f x {x'-t) 



dx' / d*v ( *-*W(» , y;t,< / )&0 B ';O 



-oo 

f 



J dr' J dt' e< ¥ -\ ee (r, r'; t, t')(3(r'; t') 
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where g2(r,r';t) is a binary quasiequilibrium distribution function (2.32), and 

<p m (x, x'; t, f ) = J dr N I n {x- t)T{t, t')I n {x'- t') Qq (x N ; t'), (3.19) 



tp n£ {x, r'; t, t') = J dr N I n {x; t)T(t, t')I e (r'; t') 6q (x N ; t'), (3.20) 
ip Sn {r, x'- 1, if) = J dr N I s {r; t)T{t, t')I n {x'; t')g q {x N ; t'), (3.21) 
m (r, r'- 1, t') = J dr N I s {r- t)T(t, t')I £ {r'- t')g q (x N ; t') (3.22) 

are generalized transport kernels which describe the kinetic and hydrodynamic dis- 
sipative processes in a system. We obtained generalized transport equations for a 
nonequilibrium one-particle distribution function and the average energy density, 
which describe both strong and weak nonequilibrium processes. From these equa- 
tions one can derive interesting limiting cases. First of all, if two last terms in the 
right-hand side of equation (|3.17| ) are not taken into account (these terms describe 
non-Markovian dissipative processes), then one obtains a generalized kinetic equa- 
tion in the mean field approximation [ || @ : 



at mar 



dx' —${\ r -r'\)—g 2 {r,r'-t)h{^t)h{x'-t). (3.23) 



From equation (|3.23|) at g2(r, r'; t) — 1 we have a usual Vlasov kinetic equation for 
a nonequilibrium one-particle distribution function. The kinetic equation ( |3.23| ) for 
fi(x; t) must be complemented by an equation for the average energy density ( |3.18|) , 
without taking into account two last terms which describe non-Markovian dissipa- 
tive processes. It is important to remember here that the binary quasiequilibrium 
distribution function g2{r, r'; t) is a functional of the local values of inversed temper- 
ature /3(r;t) and average density n(r;t): g2(r,r';t) = g2(r,r'\(3(t),n(t)). If spatial 
correlations connected with an interaction between a separated group of particles 
and the media are neglected (that is valid in the case of small densities), i.e. formally 
putting for £> q (x N ; t) in (|2.22f) the density of interaction energy £ mt (r) to be equal 
to zero (UN{r N ;t) = 0), then ( |2.22| ) transforms into a quasiequilibrium distribution 



3|) that corresponds to the usual boundary conditions of the weakening of corre- 
lations by Bogolubov [ |22|] for the solution of the Liouville equation. In this case, 
from the system of equations (|3.17| ), (|3.18|) , we have a generalized kinetic equation 



for the nonequilibrium one-particle distribution function which was obtained earlier 



by Zubarev and Novikov [ 52 



dx' [ dt' e^'^Vnn^z';*,*')^';*')- (3-24) 
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However, in this equation, contrary to the kinetic equation for /i(x; t) of paper [ |52|, 
we used a projection procedure for the exclusion of time derivatives for parameters 
a(x; t) which are determined from the self-consistency condition (|3.5| ) and in our case 



according to ( pJ2|) , ( FX22|) , (|2T23| ) and u = e, they are equal to e -a ^ = f x {x]t)/e. 



In the kinetic equation (|3.24j) , the transport kernel has the following structure: 



tp nn (x,x';t,lf) = J dr N (l-^ / (t))A 1 (x)T(t,t')(l-^'(0)^i(^)fl 



fl(Xj',t) 



3 

f 



where the projection operator &'(t) is defined as 



e 

(3.25) 



, ~ , - , * f d(A(x'))i , 
&»{t)A{x') = (A(x'))l + Jdx d yjp {h^x) - (n^Y) , (3.26) 

and the average values (A(x')) t (l are calculated with the help of the quasiequilibrium 
distribution function (|2.5| ). The kernel (p nn (x,x';t,i/) describes kinetic processes in 
the system. 

A coupled system of equations for the nonequilibrium one-particle distribution 
function and the density of the total energy (|3.17|) , ( p,18|) is strongly nonlinear. 
It takes into account complicated kinetic and hydrodynamic processes and such a 
system can be used for describing both strongly and weakly nonequilibrium states. 

In the next subsection we apply equations ( p. 171 ), (|3.18|) to the investigation of 
nonequilibrium states of the system, which are close to equilibrium. In this case the 
equations are simplified significantly. 



3.3. Kinetics and hydrodynamics of nonequilibrium state near equilibrium 

Let us assume that the average energy density (£(r))*, the nonequilibrium one- 
particle distribution function fi(x;t) and parameters P(r;t), b(x;t) deviate slightly 
from their equilibrium values. Then, being restricted to the linear approximation [ 
||, the quasiequilibrium distribution function £ q (j^ji) (|3~3|) can be expanded over 
deviations of the parameters /3(r;t), b(x;t) from their equilibrium values: 



g q (x N ;t) = Q (x N ) 



dr S/3(r;t)S(r) 



dx Sb(x; t)hi(x) 



(3.27) 



where Qo(x 

/dr, 



is an equilibrium distribution of particles, Z 



n e 



is a grand partition function, (3 = is an equilibrium value 

for inverse temperature, k-Q is the Boltzmann constant, 5/3(r;t) = /3(r;t) — (3, 
5b(x; t) = b(x; t) + /3/j,, fi is an equilibrium value of the chemical potential. In formula 
( p.27[ ) it is convenient to transform the dynamical variables S(r), fii(x) = n 1 (r,p) 
and parameters Sj3(r;t), Sb(x;t) into Fourier- components. Then we obtain: 



£ q (x N ;t) = g {x N ) 



1 - ^2 8p-h(t)£k ~/^2 d P Sb -k(p; t)n k (p) 



(3.28) 
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where £' fc = £ fc(fc ^ 0) and 



4 = / dr e- ifc r ^(r), 
™fc(p) = / dr e~ ife,T ni(r,p), 



5/3_ fe (i) = / dreP* r 6l3{r;t), 
?_ fc (p;t) = / dr e ifc r 56(x;t). 



Using the self-consistency conditions (|3.5|), let us consequently exclude parameters 
5f3_ k (t), 5b_ k (p;t) in equation ( f3.28| ). From (nfc(p))* = (n k (p)y one finds: 



dp' $ k (p,p')5b- k (p';t) = -(n fe (p))* - (n k (p)£- k )oSf3k(t) } 



(3.29) 



where (• • -)o denotes averaging over the equilibrium distribution, namely (. . .)o = 
Jdr N ... g (x N ), (n h {p)) = 0, {k ^ 0); 

= (n k (p)n-h(p'))o = n5(p - j/)/ (p') + n 2 /o(p)/o(p)^(fc) (3.30) 

is an equilibrium correlation function, n = N/V, fo(p) = (/3 /2nm) 3 / 2 e~ l3p2 ^ 2m 
Maxwell distribution, h,2(k) is a Fourier- image of the function g2{R) — 1: 



is a 



M*0 



dile-* w g 2 (fi)-l , 



(3.31) 



g2(R) — g2(\f — f'\) is a pair equilibrium distribution function. Now let us define 
the function $^ 1 (p, p') inversed to the function $ k (p, p') by means of the relation 



J dp' W^Cp'.p) = 5(p"- P ). 



(3.32) 



Taking into account ( [3.30| ), from ( |3.32|) we find $ fc 1 (p",p / ) in an explicit form 

$ k \p",p') = 5 ^';j' ) -C2(k), (3.33) 



where C2(k) denotes a direct correlation function which is connected with the corre- 
lation function h 2 {k) as: h 2 (k) = c 2 (fc)[l — ra^fc)] -1 . Let us multiply the left-hand 
side of ( |3.29|) by the function $^ 1 (p",p / ) and integrate it with respect to p'. Taking 
into account ( |3.32| ), we find: 

5b- k (p;t) = - J dp' ^(p^OMpO)*- W) / V <$> k \p,p ! ){h k (p')£_ k )«. 



Substituting the obtained value of the parameter 5b^ k (p; t) in ( 3.28|) we find after 
some transformations: 



g q (x N ;t) = Qq{x n ) 



1 - J^5^ k (t)h^ + ^ I dp dp' (n k (p'))V(p',P)n k (P 

k k J 

(3.34) 
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where 



/4 nt 



S k - dp dp' {E k h_ k (p'))^ k \p\p)h k {p) = Sf - (Sfn_ k ) S^(k)n k 



1 N 

4* = - 9 E *( 



(3.35) 
(3.36) 



N N 

|r^|)e- ifcT \ n k = J2^ ikri 

are the Fourier-components of densities for the interaction energy and the number 
of particles, respectively. Further, it is more convenient, instead of the dynamical 
variable of energy £ k , to use the variable h 1 ^ 1 (|3.35|) which is orthogonal to h k {p) by 
means of the equality: 

(/4 nt n- fe ) = 0. (3.37) 

From the structure of the dynamical variable h 1 ^ ( [3.35 ) it can be seen that it cor- 
responds to a potential part of the Fourier-component of the generalized enthalpy 
h k , which is introduced in molecular hydrodynamics [ [5], [77], [79| : 



h k = £ k - (£ k h- k ) h k = h k n + h 



k 5 



where 

h^ = £^-(£^h_ k ) S^(k)h k 
is a kinetic part of the generalized enthalpy, 



(3.38) 
(3.39) 



c k 



kin 



N 2 

i=i 



is the Fourier-component of the kinetic energy density. S2 = (n k n^ k ) is a static 
structure factor of the system. Taking into account the orthogonality of the dynam- 
ical variables h 1 ^ and n k (|3.37| ), from the self-consistency condition (/i^*)* = (h^Ya 



that is equivalent to (£*•)* = (£k)\ one defines the parameter 5f3 k {t): 

8W) = -C^y^Kk), 



(3.40) 



where $ hh (k) = (h^ti^Q. Finally, let us substitute (|3.40|) into ( |3.34|) . As a result, 
we obtain a quasiequilibrium distribution function in the linear approximation: 



Qq {x ;t) = Q (X ) 



1 + E'^v^)^ + E7 dp dp ' (^(P'))W,P)n k (p) 

k k J 

(3.4 



Taking into account (|3.41|) , the nonequilibrium distribution function g(x N ;t) ( |3.10|) 

(3.42) 



or (|3.12|) has the following form in this approximation 

g [x N ; t) = Qo{x N ) x 

"i+£'#?>**S(*)£S t + 
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E' I dp dp' (n k (p')}V(P,P)Mp) ~ 

k 

t 

E' / ^^-^{^^(fc^^oirc-fc)- 

K -oo 

t 

E' / d ^ d P' / dt'e^-*)(n fe (p'))^fe 1 (p / ,p)ro(t,i , )4(-fe,p) 



where 



/„(*!, p) = (1 - ^oK(p), 
if(fc) = (1 - ^Aj* 



(3.43) 
(3.44) 



are generalized fluxes in the linear approximation, n k (p) = iL^n k {p), h 1 ^ = iL^h 1 ^, 
To(t,t') = e( t ~ t ')( 1 ~^' ) lLjv is a time evolution operator with the projection operator 
which is a linear approximation of the generalized Mori projection operator 
3P{t) ( |3.15|) , with taking into account orthogonalization of variables {S k ,n k (p)} — > 
{h 1 ^, n k (p)}. According to the structure ( |3.41| ), ^ acts on the dynamical variables 

0> O A W = ^(A«h™l) ^(k)ti* + E / ^ dp' (i fc m_ fc (p')) $ fc 1 (p',p)n fc (p). 



(3-45) 

For kinetics and hydrodynamics of nonequilibrium processes which are close to 
an equilibrium state, the generalized transport equations ( |3.17| ), ( |3.18| ) in the linear 
approximation for the nonequilibrium distribution function g (x N ;t) ( |3.42 ) trans- 
form into a transport equation for f k (p;t) = (Wfc(p))*, h^it) = {h 1 ^) 1 

«i/*(p;t) + — /fc(P5*) = -— -nf (p)c 2 (k) / dp' f k (p';t)+in nh (k,p)^(t)- 
at m m J 

t t 
J dp' J dt' e e (*'- t Vnn(fe,P,p';t,0/fc(P , ;0- / M e s V-»<p nh (k,p;t,t')h}*(t'), 

— oo — oo 

(3.46) 



d_ 
dt' 



h?(t) = / dpin hn (k,p)h(p;t)- 



(3.47) 



J dp J dt' e^'' t ^ hn (k,p;t,t')f k (p;t')~ J dt' e E <*-*W(M,0>ik t (0, 

— oo — oo 

where iQ n h(k,p), ifi/ in (fe,p) are normalized static correlation functions: 

m nh (k, P ) = (h h (p)hp h ) $£( k ), (3.48) 



: int 



ifl hn {k,p) = / dp' (h k n_ k {p'))^ k \p',p 



(3.49) 
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and 

<p nn {k,p,p';t,t') = I dp" (I n (k,p)T (t,t')I n (-k,p")) ^ k \p",p'), (3.50) 

Vhn {k,p-t,t') = J dp' (If(k)T (t,t')I n (-k,p')) ^ k \p',p), (3.51) 

Vnh {k,p-t,t') = (I n (k,p)T (t,t')ff{-k)) <f>^(k), (3.52) 
Vhh (k;t,t') = (l£\k)T (t,t')If(-k)) $- h l(k) (3.53) 

are generalized transport kernels (memory functions) which describe kinetic and 
hydrodynamic processes. The system of transport equations (|3.46|) , ( |3.47|) is closed. 



Eliminating h^(t), it is possible to obtain a closed kinetic equation for a Fourier 
component of the nonequilibrium one-particle distribution function. For this purpose 
one uses a Laplace transform with respect to time, assuming that at t > the 
quantities f k (p] t — 0), h^(t = 0) are known 



A(z)=i dte izt A(t), z = u + ie, e -> +0. (3.54) 
Jo 

Then, equations (6.46) and (6.47) are presented in the form: 

zf k (p; z) + ^/ fc (p; z) = -^n/ (p)c 2 (fc) / dp' Up'; z) + (3.55) 
m m J 

E nh (k,p; z)h k \z) - J dp' ip nn (k,p,p'; z)f k (p'; z) + / fe (p; t = 0), 
zh™(z) = J dp' E hn (k,p'; z)hip'] z) - m (fe; z)K\z) + hf(t = 0), (3.56) 



where 



E nh (k,p; z) = iQ nh (k,p) - (p nh (k,p; z), (3.57) 
T, hn (k,p; z) = iQ hn (k,p) - (p hn (k,p; z). (3.58) 



Let us solve equation ( |3.56| ) with respect to h^i^z) and substitute the result into 



( P . 5 5|) . Then one obtains a closed kinetic equation for f k (p; z) (at h x k {t = 0) = 0): 

\k ' p lk m p f 

zfk(p; z) H f k (p; z) = nf (p)c 2 (k) / dp' f k (p'; z) - 

m m J 

J dp' D nn (k,p,p';z)f k (p';z) + f k (p;t = 0), (3.59) 

where 

D nn (k,p,p';z) = ip nn (k,p,p';z) - T, nh (k,p; z) — ■ ^— — rE hn (k,p'; z) (3.60) 

z + (phh{k;z) 

is a generalized transport kernel of kinetic processes, which is renormalized taking 
into account the processes of transport of the potential energy of interaction between 
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the particles. If we put formally in D nn (k,p,p'; z) that h mt k = (such equality is 
valid when the contribution of the average potential energy of interaction is much 
smaller than that of the average kinetic energy), then from ( ft. 59 ) one obtains a 
kinetic equation for /fc(p; z), 



zfk(p; z) + — ^/ fc (p; z) 



ik ■ p 



m 



in 



nf (p)c 2 (k) / dp' f k (p';z 



J dp' <p' nn (k,P,p'; z)h( P '; z) + f k (p; t = o). 



(3.61) 



Equation (|3.61|) was derived for the first time by using the Mori projection operators 
method in [ |107] , |108| , |109|| , when the microscopic phase density n k {p) was a param- 
eter of an abbreviated description. In this case the memory function ip' nn (k,p,p'; z) 
has the following structure: 



where 



fc,p,p';M') = / dp" (I n (k,p)T^t')I° n (-k,p")) ^ 



J°(fc,p) = (l-^K(p) 



k' 



P",P' 



(3.62) 



(3.63) 



is a generalized flow, is the Mori projection operator introduced in [ |107| , 110 



1091 : 



^' Ay = J2f d P d P (A k ^ k (p'))o^ k 1 (p',p)h k (p), (3-64) 

k 

Tq(£, t') = e^~*'^ 1_ ^'o) lL A r is the corresponding time evolution operator. The kinetic 
equation (|3.59|) for the Fourier component of the nonequilibrium one-particle distri- 
bution function stimulates investigations of the dynamical structure factor S(k;ou), 
time correlation functions for transverse and longitudinal flows of the particles, dif- 
fusion coefficients, viscosities for dense gases and liquids [ [77], [TJ| [TDTJ |108| , |109| , 
TO [TTT1 , Imj [TT3| , [TT|, [TT5|, EH, Pl7l, PI. In papers by Mazenko, there was 
performed a derivation of the linearized Boltzmann-Enskog equation by expanding 
memory functions ip' nn (k,p,p';t,t') on density. In the case of weak interactions a 
Fokker-Planck equation was obtained. Mazenko [ |114| and Forster [ |115| , |136|1 cal- 
culated the generalized viscosity coefficient with the help of kinetic equation ( |3.59 ). 
Self-diffusion of a one-component plasma [ |181|| , as well as spectra of mass and charge 



fluctuations in ionic solutions were investigated on the basis of this equation. In ac- 
tual calculations, certain approximations for the memory functions ip' nn (k,p,p'; t, t') 
were used. However, the main drawback of the kinetic equation (|3.59|) lies in its non- 



self-consistency with the conservation law of the total energy for dense gases and 
liquids when the contribution of the interaction energy between the particles to 
thermodynamic quantities and transport coefficients becomes important. This fact 
was also pointed out earlier in the paper by John and Forster [ |118[ , where they 



performed an investigation of the dynamical structure factor S(k; oj) in the inter- 
mediate region of wavevector k and frequency values u for a simple liquid on the 
basis of the set of dynamical variables of an abbreviated description n k (p) and E k . 
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In the next subsection, on the basis of a system of transport equations for Fourier 
components of the nonequilibrium one-particle distribution function and the poten- 
tial part of enthalpy ( [3.46 ), Q3.47Q , we shall obtain equations for time correlation 
functions. We shall also investigate the spectrum of collective excitations and the 
structure of generalized transport coefficients. 

3.4. Time correlation functions, collective modes and generalized transport 
coefficients 



With the help of combined equations (|3.46|) , (|3.47| ) one obtains a system for time 
correlation functions: 



$ nn (k,p,p';t) = J dp" (h k (p;t)n_ k (p";0)) ^ k 1 (p",p') 

® hn {k, P] t) = J dp' (^ t (t)n_ fc (p / ;0)) $ fc 1 (p / ,p), 

$ nh (k,p;t) = (n fe (p;t)^(0)}o^ 1 (fe) ) 
«Mfc;f) = K\t)ti«i(0)) o ^(k), 

where h k (p;t) = e-' lLNt n k (p; 0), h^(t) = e~ iIwt /ijf (0). 
One uses the Fourier transform with respect to time 



(3.65) 

(3.66) 

(3.67) 
(3.68) 



(a) w = / dte^(ay. 
Then we write the system of equations ( |3.46| ), ( |3.47| ) in the form: 
- iu(h k (p)) w = dp' E nn (k,p,p';uj + ie)(n k (p')) UJ + E nh (k,p;u + ie)(/4 n %,(3.69) 



iw(h^) u = dp' E hn (k,p';u + ie)(h h {p'))„ - <p hh (k;u + ie){h}^) UJ , (3.70) 



(3.71) 
(3.72) 
(3.73) 



where 

S„ n (fc,p,p';u; + ie) = ifi nn (fe, p, p') - ip nn {k, P, p'\ w + ie) 
J: nh (k,p;uj + ie) = iQ nh (k,p) - cp nh (k,p; u + ie), 
E hn (k,p;uj + ie) = iQ hn (k,p) - (p hn (k,p; u + ie). 

It is more convenient to present the system of equations ( |3.69| ), ( |3.70j ) in a matrix 
form: 

- iu;(a fc ) w = E(fc; u + ie) (dk)^, (3.74) 
where d k = col(n k (p), h 1 ^) is a vector-column and 

/ dp' T, nn (k,p,p'; u + ie) E nh (k, p, u + ie) 
j dp' E hn (k,p';uj + ie) -ip hh (k; uj + ie) 



S(fc;cu + ie) 
t(k;uj + ie) 



(3.75) 



dt e^ +i£)t t(k;t). 
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Now, one uses the solution to the Liouville equation in approximation (|3.41|) without 
introducing the projection operator & q {t): 

t 

g (x N - 1) = g q (x N ; t)- j dt' tfV-WvV-Q + iL ^j g q (x N ; t>) . 

— oo 

Then, from the self-consistency conditions (afc)* = (ak) q one obtains a system of 

equations which connects the average values {hk(p))u and (h 1 ^)^ with spectral func- 
tions of time correlation functions: 



i(uj + ie) 



J dp' <£> nn (k,p } p';uj + ie) $|$(fc,p; u + ie) 
j dp' <(fc,p'; cu + ie) u + ie) 

J dp' $ nn (fe,p,p') 
-$hh(k) 

jdp' <$> nn (k,p,p';uj + ie) ^(k,p;u + ie) 
J dp' uj + ie) -n£(k; u + ie) 





(n k {p)}u 


X 







(n k {p))u 


X 





Such a system can be presented in the next compact form: 

ia;$(fc;u; + ie)(a fc ) w = $(fe) - i(w + ie)<§(fc; w + ie) (a fc ) w . (3.76) 



Let us multiply equation ( 3.74TJ by the matrix $(fc; co> + ie) and compare the result 

(3.77) 



with equation (|3.76|) . So we find 

2$(fc; 2) = E(jfe; z)<£(fc; z) - <f>(fc) 

/■oo 

$(fc;z) = i / dt e izt $(k;t), 



z = uj + ie, 
e^O 

Jo 

or in an explicit form: 

z$ nn (k,p,p';z) = J dp" Z nn (k,p,p";z)<f> nn (k,p",p';z) + (3.78) 
E nh (k, p, (fc, p'; z) - $ nn (fe, p, p')> 

dp ,/ s nn (fe J p J p , V)<t(fe,P /, ^) + Sn /l (fe,P^)^(fc;^), (3.79) 

dp" S ftn (fc,p";z)^(fc,p^p^) - m (fc;z)^(fc,p';z),(3.80) 

dp" S hn (fc,p";z)^(fc,p"^) - ip hh (k-z)^(k-z) - $ hh (k), 

(3.81) 

where the condition Qh n (k,p') = § n h(k,p) = is taken into account. 



z$^(k,p;z) 
z^(k,p';z) 
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Therefore, we obtained a system of equations for normalized correlation func- 
tions ( |3.65| )-( |3.68| ) for the description of nonequilibrium states of the system when 
kinetic and hydrodynamic processes are considered simultaneously A similar system 
of equations for time correlation functions in the method of projection operators was 
obtained in the paper by John and Forster [ |118|| . In this paper the parameters of 
an abbreviated description were orthogonal dynamical variables hk{p) and So- 
lutions of such an equation system in the hydrodynamic limit were found using the 
projection procedure for the first five moments of the nonequilibrium one-particle 
distribution function. In view of this, the dynamical structure factor S(k\ u) for 
intermediate values of k and u was investigated using a parametric approximation 
of memory functions for liquid Argon. However, in this paper the question con- 
cerning collective modes and generalized transport coefficients, when kinetics and 
hydrodynamics are connected between themselves, is not discussed. 

In order to solve the system of equations (|3.78|) — (|3.81 ) we also apply the pro- 



jection procedure [ 11081. Let us introduce the dimensionless momentum £ = — 



moo ' 



Vq = (m/3) 1 . Then the system of equations ( |3.77j ) can be rewritten in the matrix 
form: 

z$(k- i, z) - E(fc; £, z)$(fc; £" , z) = -$(fc; & £'), (3-82) 

where it is clear that the integration must be performed with respect to the repeating 
indices Further, let us introduce the scalar product of two functions, and 
V>(£), as 

<M) = j d£<f(£)/o(£M£)- (3-83) 
Then, the matrix element for some "operator" M can be determined as 

(0|M|V> = j d£ d£' ^«)M(£,0/o(0^(0- (3-84) 

Let <f>(£) = {0,u(£)} be the orthogonalized basis of functions with the weight /o(£), 
so that the following condition is satisfied: 

(0,1^) = ^, ^2\(f>u)(<f>u\ = l, (3.85) 

V 

where 

MZ) = <t>i™n(€) = (llmlnl)- 1 / 2 H t (Q H m (Q H n (Q , (3.86) 

Hi(£) = 2~ l / 2 Hi(£ i /2), Hi{£) is a Hermite polynomial. Then, each function in the 
system of equations (|3.82|) , which depends on momentum variables £, can be 
expanded over functions in the series: 

z) = X>;(*)s*m(M^(0/o(0, (3-87) 

nk-^i'-z) = J2^)K(k;z)MOfo(C), (3.88) 
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where 

$^(fe;z) = (^|$(fe;^^z)|^) =J d^d^^/oCO^fc;^^^)^). (3-89) 
£„„(M = <^|£(fc;^';z)|</v> =/ d£ d£' #(*)/o(OE(fc; & ( 3 - 90 ) 



Let us substitute expansions (|3.87| )- (|3.90| ) into equation (|3.82|) . As a result, one 
obtains: 

z$„ M (fc; - X) *)$7/*( fe ; *) = -*M*0- ( 3 - 91 ) 

7 

In actual calculations, a finite number of functions from the set <fr u (£) is used. Taking 
into account this fact, let us introduce the projection operator P which projects 
arbitrary functions onto a finite set of functions 



(3.92) 



u=l 



v=l 



Here n denotes a finite number of functions. Then, from (|3.91|) we obtain a system 
of equations for a finite set of functions At (^), 



^ U^ 7 - ifVy(fc) + A^(fc; 

7=1 



$ 7M (fc;z) = 



(3.93) 



where 



■D Vfl (k] z 



0{k;z) + E{k;z)Q z-QE{k;z)Q QE{k; (3.94) 



are generalized hydrodynamic transport kernels and 

i^(fc) = (^|ifi(fc)|0 M ) 



(3.95) 



is a frequency matrix. Note that matrices ifl(fe) and <^(fc; z) are defined according 
to (|3~48D , (Ertgp and (^50|)-(|3~5l). 

Let us find solutions to the system of equations (|3.93|) in the hydrodynamic region 
when a set of functions M (£) present five moments of a one-particle distribution 
function: 



1 



0i(O = i, 2 (o = e 

Then, the following relations are fulfilled: 

(l|n fe (£) = / d£ n k (g) = n k 



-,{£) = $,. (3.96) 



(kM*) = / d£Mt)$r=& ( 3 - 97 ) 

(6-1/2(^2 _ 3)]M ^ = y ^ . fc(£) g-1/2^2 _ 3) = ^kin _ 3^-1 = ^to 
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for the Fourier components of densities for the number of particles, momentum and 
the kinetic part of generalized enthalpy. Besides that, the microscopic conservation 
laws for densities of the number of particles and momentum can be written in the 
form: 

(l|"ft(£) = -ikyPim- 1 , 

- (3.98) 



<fa|Wfe(£) 



where TU a is a Fourier component of the stress-tensor. 

If we choose the direction of wavevector k along oz-axis, then </>„(£), v = 1, 2, 3 
will correspond to longitudinal modes, whereas </>„(£) at v = 4, 5 will be related to 
transverse modes. 

From the system of equations fl3.93Q , at v — 4,5, 04 (^) = 0s(£) = £ y , one 
obtains an equation for the Fourier component of the time correlation function 
connected with the transverse component of the momentum density $44(^5 z). From 
this equation one finds: 



$ 44 (A;; z) = $fjfc; z 



1 



where 



-jjp 

^i (kin) (fe;^) 



D^)(fc;*)+i^>(fc;*)> 



E(fe, *)Q z - QS(fe, z)Q QE(fc, £, z) 



-1 



D (k;z) = ik 2 r](k; z)(mn) 1 , 



(3.99) 



(3.100) 
(3.101) 
(3.102) 

nn 

(3.103) 
(3.104) 



where 77 (fe; z) denotes a generalized coefficient of shear viscosity. Such a coefficient 
consists of two main contributions. The first one is D^p (fe; z), whereas the second 
contribution Dpp int \k; z) describes a relation of kinetic and hydrodynamic processes. 
If the last term is neglected, which formally corresponds to h 1 ^ = 0, then one obtains 
an expression for the generalized coefficient of shear viscosity 77 (fe; z), obtained earlier 



by the authors of [ 114, |115|| when solving the equation 



z$ nn {k,p,p';z) - j dp" iQ nn (k,p,p")$ nn (k,p",p';z) + 



(3.105) 



j dp" ip' nn {k,p,p";z)$ nn {k,p",p';z) = -<$> nn {k,p,p') 



by using the method of projections [ |108|| in the hydrodynamic region. Equation 
( |3.105| ) corresponds to kinetic equation ( 3.61 ) for the Fourier component of the 
nonequilibrium one-particle distribution function. In (p' nn (k,p,p'; z) the projection 
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is performed on a space of dynamics of the "slow" variable hk(p) connected with 
the microscopic conservation laws for the particle-number density and momentum 
( p.98|) . In our approach in Dp p (k;z), the projection is carried out on an extended 

space of dynamics of the "slow" variables n^p) and h 1 ^. Moreover, h 1 ^ is a hydro- 
dynamic variable, whereas n^p) is a kinetic one and only its average value with five 
moments ( |3.96| ) can be related to hydrodynamic quantities. 

If we put v = 1, 2, 3, <M£) = 1, M£) = Mi) ~ Q- l! \e ~ 3) in the system 
of equation ( |3.93|) , then we obtain: 



+ 



z® na (k; z) - iQ np (k)$ pa (k; z 

z% a (k; z) - ia pn (k)<5> na (k; z) + L>JJ p (fc; z)% a {k; z 
E ph Mn(k; z)$ fc u» (fc; z) - S pftint (fc; z)$ fcta t a (fc; z 

z® h Kn a {k; z) - E h un p (k; z)<S> pa (k; z 
D h km h *in(k; z)® h Mn a (k; z) + D h ^ K ^{k\ z)$ h mt a (k; z 

z® hi nt a (k; z) - E fttatp (fe; z)% a (k; z 
D hi nt h ^n(k; z)$ h un a (k; z) + D h i^ h int(k; z)^^nt a (k; z 

where a = {n fc ,p fc , /i| in , h 1 ^} and 



E p ^kin(fc; z) — iQ ph kin(k) — D ph kin(k; z), 

E pfe tat(fe; z) = iQ ph int(k) - D ph int(k; z), 

E ft km p (fc;2) = ii\km p (fc) — D h kin p (k] z), 

Y, h im p (k; z) = iQ h int p (k) - D h int p (k; z), 



-$ na (fc), (3.106) 
(3.107) 

(3.108) 

-$ ft kin a (fc), 

(3.109) 

-® h M a (k), 



(3.110) 



iQ a b(k) and D a b(k; z) are determined according to (|3.95| ), ( p.94|) . From the system 
of equations (|3.106|) -( |3.109|) one can define the Fourier components of the particle 
number density correlation functions 

$ nn (k;z) = $n(fc;2) = (l|$ nn (fc, £, , (3.111) 
as well as of the longitudinal component of the momentum density 

$ p [ p (k;z) = $ 22 (fc;z) = (t z \®nn(k,Z,£';z)\Q, (3.112) 

for the kinetic part of generalized enthalpy 

c^Mn^fc; z) = $ 33 (fc; z) = (Q- l/2 (e - 3)|$ nn (fe, £, ^-^((O 2 - 3)) (3.113) 

as well as for the potential part of generalized enthalpy $^int^mt(fe; z) and cross 
correlation functions, especially ^^^(fe; z), $ n ^kia(fe; 2), $ n/l mt(fc; 2), $ p/i km(fc; 2), 
$p^tat(fe; 2). It is important to point out that the system of equations ( |3.106| )- (|3.109|) 
corresponds to the system of equations for Fourier components of the average values 
of densities for the number of particles (rifc) 2 , longitudinal momentum (Pk) z , kinetic 
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(h^ n ) z and potential {h^) z parts of generalized enthalpy: 



z(h k ) z - itt np (k)(p k ) z = 

z(p k ) z - \Vt pn (k){n k ) z + 
Dl(k; z)(p k ) z - S p , kin (fc; z)(h^) z - E ph ^(k; z)(h^) z = 

z(tit) z -Y; h u« p {k-z)(p k ) z + 
IWd»(fc; z)(h^) z + D h u^(k; z)(h^) z = 
z{U^) z -E h i a t p {k;z)(p k ) z + 
iW^fc; z){tit) z + £W*(fc; z)(h^) z = 



-(n k (t = 0)), (3.114) 
(3.115) 

-<p*(* = o)>, 

(3.116) 

-(^ in (t = o)), 

(3.117) 



This system of equations is similar in construction to the equations of molecular hy- 
drodynamics [ [nj. The difference consists in the fact that instead of the equations 
for the Fourier component of the mean enthalpy density (h k ) 1 which is introduced in 
molecular hydrodynamics, there are two connected equations for the Fourier compo- 
nents of mean values of the kinetic and potential parts of enthalpy density. Moreover, 
instead of the generalized thermal conductivity which appears in molecular hydro- 
dynamics, the dissipation of energy flows is described in equations (|3.114| )- (|3.117|) 
by a set of generalized transport kernels Z\km fe km(A;; z), D /l ki nh int(fc; z), D^m^i^k] z), 
D h int h int(k; z). Obviously, transport kernels give more detailed information on the 
dissipation of energy flows in the system because they describe the time evolution of 
dynamical correlations between the kinetic and potential flows of enthalpy density. 

Solving the system of equation ( |3.106| )- (|3.109|) at a = n, one obtains an expres- 
sion for the correlation function "density-density" $ nn (fc; z) 



<& nn {k-z) = -S 2 (k) 



i^np(fe)i^pn(fc) 

z + Dpp(k; z) 



-i 



(3.118) 



where 



D n>( k '> z ) = D p P ( k 'i z ) ~ Z ph *in(k; z) [z + D h un h un(k; z)] 1 E ft ,kin p (fe; z) - 

E ph mt(k;z) [z + D h int h mt(k; z)]' 1 T; h m p (k; z), (3.119) 

E ph un(k;z) = E ph u*(k; z) - (3.120) 
S p/l int(fe; z) [z + D h int h int(k; z)}~ D h int h Mn(k; z), 

(3.121) 



S^kin^-, ( k , Zj — /,kin T) ( k , Z 



/l ki 



l pV 1 



D h um h int(k; z) [z + D h int h mt(k; z)] 1 z) 



D h kin h kin(k; z) = D h kin h kin(k; z) — (3.122) 

D h kin h int(k] Z) [Z + D h i-at h ^{k] Z)}~ 1 D h int h kin(k] z). 



In expressions (|3.119|) - (|3.122|) we can observe an interesting renormalization of the 
functions £ a b and via the generalized transport kernels for fluctuations of flows 
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of the potential part of enthalpy density. It is important to point out that in the 
mode-coupling theory developed by Gotze [ |95], |96|, an expression for $„„(&; z) has 



the same form as in ( FTTTgp . However, D p l p (k; z) is connected only with the gen- 
eralized shear viscosity r)"(k;z), since the densities of the number of particles 
and momentum p k are included in the set of variables of an abbreviated descrip- 
tion. In our case D pp (k; z) takes into account both thermal and viscous dynamical 
correlation processes. Excluding from (|3.118| ) the imaginary part <&|| n (fe; a>) of the 
correlation function $ nn (fe; z), one obtains an expression for the dynamical structure 
factor S(k; uj) in which contributions of transport kernels corresponding to the ki- 
netic and potential parts of the enthalpy density are separated. It is evident that 
the main contribution of the transport kernel D h mt h mt(k; z) to the S(k; uj) for liquids 
was in the hydrodynamical region (the region of small values of wavevector k and 
frequency uj), whereas D^^ik] z) will contribute to the kinetic region (orders of 
interatomic distances, small correlation times). In the region of intermediate values 
of wavevector k and frequency u, it is necessary to take into account all the transport 
kernels S p ^km(fc; z), Splint (k; z) , D h kin h kin(k; z), D h int h kin(k; z), D h int h mt (k; z) . Since it 
is impossible to perform exact calculations of the described above functions, it is 
necessary in each separate region to accept models corresponding to the physical 
processes. Obviously, it is necessary to provide modelling on the level of general- 
i zed t ra nspo rt kernels (p nn (k,p,p r ;t,f), <p nh (k,p;t,t'), <p hn (k,p r ;t,f), <p hh (k;t,t r ) 
( p.50| )-( p.53| ). In particular, we can accept for <p nn (k,p,p';t,t') the Boltzmann- 
Enskog model for hard spheres or the Fokker-Planck model for the case of weak 



correlation, as in the papers by Mazenko [ 1 11, 114 , 135 1. Such approximations can 



be applicable to kinetic processes at interatomic distances (region of large values 
of wavevector k) and small correlation times (region of large frequencies uj) of in- 
teratomic collisions. Besides, it is necessary to take into account information about 
interactions of particles at short distances. In the transport kernels Lp nn {k; t, t'), 
<finh{k; t, f), iphn(k; t, t'){phh(k; t, t') it is important to distinguish processes con- 
nected with a type of particles interaction at short and long distances, analogously 
to the case of Enskog-Landau kinetic equation for a model of charged hard spheres. 
The problems of modelling for these kernels are complicated by the absence of a 
small parameter in terms of which the perturbation theory might be developed. 



The transport kernel ip' nn (k, p, p'; t, t') ( |3.62|) and its moments on space [l,pi, f-] 
were analyzed in the limit z — > oo and the hydrodynamic limit k — > 0, z — > 



in the paper by Forster [ LIS]. The modelling problems of transport kernels for 
intermediate values of k and u are reflected in the details of the description of 
spectra for collective excitations, as well as in the dynamical structure factor. In 
the limit k — > 0, uj — > 0, the cross correlations between the kinetic and potential 
flows of energy and shear flows become not so important and the system of equations 
( 3. 1 IS )— ( 3.122| ) gives a spectrum of the collective modes of molecular hydrodynamics 
77], [75[ . For intermediate values of k and u, the spectrum of collective modes can 
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be found from the condition 



z if2 n p(fc) 
i£l pn (k) z + -Dp P (fc; z) 

Zjfjkaxp 
S/jintp 



X] h int ( & ^ ZJ 





k; z) z + D h kin h kin(k; z) D h kin h int(k; z) 
D h int h kin(k; z) z + D h int h int(k; z) 



S p/l kin(fc; z) 



0, (3.123) 



in which contributions of the kinetic and potential parts of generalized enthalpy 
are separated. This will be reflected in the structure of a heat mode at concrete 
model calculations of the wavevector- and frequency-dependent transport kernels 
D h kin h kin(k; z), D h kin h int(k; z), D h int h kin ( k ; z) , D^iat^int (fe ; z) depending on k and to. 

The system of equations ( p . 9 3|) for time correlation functions allows an extended 
description of collective modes in liquids, which includes both hydrodynamic and 
kinetic processes. Including on the basis of functions <f> v (£) (|3.96|) some additional 
functions, 

vi(*) = - 5)6, vfite) = ^-m - (3.124) 

which corresponds to a 13-moments approximation of Grad, one obtains a system of 
equations for time correlation functions of an extended set of hydrodynamic variables 

{n k , p k , h**, U k , Q k , h^y. 



where 



§ K (k;z) 



-ifi H 


(fc)l> H (fc 


;^) + ^ H 


(fc;z 


)$ H (fc;z; 


= -$H 




^np 


$„ A kin 


$„n 








^pp 




$ P n 






$^kir 


n &h kin p 


3>/iidn/ikta 




n ^/i kiii o 








$ nh kin 


$nn 




















n ^h lnt p 






n 3>h int Q 





(3.125) 



(3.126) 



(*;*) 



is a matrix of Laplace images of the time correlation functions, IIfc= fd-€ ( Pn(£)'nh(£), 



iQ u (k) 






iQ n p 














if2p n 





iJlp/jkm 


if2 p n 





in 





ifi/jkinp 










> o 





if2n> 








i^nQ 


















in 













iO/jintQ 






ph h 



(3.127) 



(fc) 
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is an extended hydrodynamic frequency matrix, 



(p a (k-z) 




D 
D 
D 
D 
D 

















H 

PP 


D H 

^ph kin 


U pTl 


u pQ 


D H 


h kin p 


-D/jkm/jkin 


D h kin U 




-C/i kin h int 


Up 


D nh kin 


Dun 


Dug 




Qp 






Dqq 




h int p 


D/jint/jkin 






D h int h int 



(3.128) 



(k;z) 



is a matrix of generalized memory functions, elements of which are transport ker- 
nels (3.94) projected on the basis of functions <j) v {£) ( |3.96|) , ( [3.124| ). For such a 
description, the spectrum of generalized collective modes of the system is deter- 



zl-itt u (k) + fi u (k) 







mined for intermediate k and to by the relation det 

which takes into account kinetic and hydrodynamic processes. In the hydrodynamic 
limit k — > 0, lo — > 0, when the contribution of cross dissipative correlations between 
the kinetic and hydrodynamic processes practically vanishes, the system of equa- 
tions for the time correlation function (|3.125|) , after some transformations, can be 
reduced to a system of equations for time corelation functions of densities for the 
number of particles n^, momentum p k , total enthalpy h^, the generalized stress ten- 
sor ii k = (1 — ^ H )iLjvPfe and the generalized enthalpy flow q k = (1 — <^H)i£jv^fc, 
where is the Mori operator constructed on the dynamical variables {n*,, p k , hk\. 
For such a system of equations, the spectrum of collective excitations is determined 
from [Pgl TO 



z 


\Vl n p 











i£lp n 


z 




1^P7T 








iQ hp 


z 











1^7rp 





z + ip n7T iQ nQ 










itt c 







0. 



(3.129) 



(*>;*) 



In the hydrodynamic limit this gives: the heat mode 

z H (k) = D T k 2 + 0(k 4 ), 
two complex conjugated sound modes 

z±(k) = ±ia;s(fc) + z s (k), 



(3.130) 



(3.131) 



where a>g(fc) = ck + o(k 3 ) is a frequency of sound propagation, Zs(k) = Tk 2 + 0{k A ) 
is a frequency of sound damping with the damping coefficient T; two nonvanishing 
in the limit k — > kinetic modes 

(3.132) 



z n (k) 
z Q {k) 



^(0) + 0(k 2 
^ QQ (0) + 0(k 2 ) 



Here Dt denotes a thermal diffusion coefficient 



Dn 



J TQ 



A 



7<^qq(°) mncp 



J TQ 



QQ 



h 2 



nc v 



7 = Cp/Cy 
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c p and cy are, correspondingly, thermodynamic values of specific heats at the con- 
stant pressure and volume, A is a thermal conductivity coefficient, h denotes a ther- 
modynamic value of enthalpy, c = ■y/mnS(k = 0) denotes an adiabatic sound veloc- 
ity. In (ETT5ID 

r=~( T - i)d t + ^ 



(3.133) 



is a sound damping coefficient with 



77 1 



p-TT 



^7T7r(0) 



4 

3 ' 



/nm, 



mg(0)$^(0)-7 
mnS(0) 



where 77 and k are shear and bulk viscosities. This spectrum coincides with the 
results of papers [ |139| , 141 , 142 1. However, at fixed values of k and ou, the trans- 



port kernels <p ww , <p W Q, Pq-k, Pqq are expressed via the generalized transport kernels 
D ull (k;z) of matrix (P^), i.e. via (p nn (k,p,p';t,t'), ip nh (k,p;t,t'), (p hn (k,p';t,t'), 



iphh(k;t,t') (|3.50|) - (|3.53D , according to the definition D vp {k\z) (|3.94|) . Here, it is 
important to point out that passing from the system of transport equations of a 
self-consistent description of kinetics and hydrodynamics to the equations of gener- 
alized hydrodynamics, we can connect the generalized transport kernels (|3.50|) - (|3.53|) 
with the hydrodynamic transport kernels D VVL (k\ z) in fl3.123j ) or fl3.128|) . Therefore, 
we gain the aim of the present section, namely, to connect transport kernels of a 
self-consistent description of kinetics and hydrodynamics with hydrodynamic trans- 
port kernels for time correlation functions for densities of the number of particles 
hk, momentum p k , enthalpy h^, which can be calculated by using the MD method [ 

1421 , HH HM Ell' B B B B B HH E03> HH> as wel1 as b y experiments 
on light scattering for different real systems. Evidently, the most important problem 
is model calculations of transport kernels (|3.50|) - (|3.53|) . Such calculations must take 
into account a type of interaction of particles at short and long distances, as well 
as their structural distribution. From this point of view, of significant interest is the 
modelling of transport kernels ( |3.50| )-( |3.53| ) for a system of charged hard spheres, 
since for a system of hard spheres a transport kernel of the type (p nn (k;p,p']t,t') 
( p.50| ) was computed via the Enskog collision integral 



120, 121 



4. Conclusion 



We have presented one of the approaches for a self-consistent description of ki- 
netic and hydrodynamic processes in systems of interacting particles, formulated 
by the NSO method. It is based on a modification of the boundary condition to 
the Liouville equation, which takes into account a nonequilibrity of the one-particle 
distribution function, as well as the local conservation laws which constitute a ba- 
sis for the hydrodynamic description. Using such a description, generalized kinetic 
equations for dense gases and liquids can be derived. The result obtained can be 
extended to quantum systems of interacting particles [ [17], [L8| . In particular, it 



can be applied to nuclear matter and chemically reacting systems, where kinetic pro- 
cesses play an important role together with hydrodynamic ones. The method applied 
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can be extended to the investigations of nonlinear and hydrodynamic fluctuations 
in polar, ionic and magnetic liquids and electrolyte solutions. 

There is a significant interest in investigations of collective excitations, time cor- 
relation functions and transport coefficients on the basis of equations ( |3.93| ), Q3.94D , 
( [3.123| ), ( 3.128|) at the presence of transport kernels tp(k; z) which describe an in- 
terference between kinetic and hydrodynamic processes. In particular, the transport 
kernel ip(k, p, p'\ t) can take into account processes of interaction of particles at short 
and long distances (at Coulomb and dipole interactions of particles). To do this, we 
can apply an approach of collective modes [ [j^, [Ll3|, [L4|, [L45|, [Ll6|, [b|7|, fb|| |TI9| , 
1501 EH S 
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